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We report on measurements of hadronic and leptonic cross sections and Icptonic forward- 
backward asymmetries performed with the L3 detector in the years 1993 — 95. A total luminosity 
of 103 pb"-*^ was collected at centre-of-mass energies ^/s « mz and ^/s « mz ± 1-8 GeV which 
corresponds to 2.5 million hadronic and 245 thousand leptonic events selected. These data lead to 
a significantly improved determination of Z parameters. From the total cross sections, combined 
with our measurements in 1990 — 92, we obtain the final results: 



An invisible width of Finv = 499.1 ± 2.9 MeV is derived which in the Standard Model yields for 
the number of light neutrino species N^, = 2.978 ± 0.014. 

Adding our results on the leptonic forward-backward asymmetries and the tau polarisation, 
the effective vector and axial- vector coupling constants of the neutral weak current to charged 
Icptons are determined to be = -0.0397 ± 0.0017 and = -0.50153 ± 0.00053. Including 
our measurements of the Z ^ bb forward-backward and quark charge asymmetries a value for the 
effective electroweak mixing angle of sin^^w = 0.23093 ± 0.00066 is derived. 

All these measurements are in good agreement with the Standard Model of electroweak inter- 
actions. Using all our measurements of electroweak observables an upper limit on the mass of the 
Standard Model Higgs boson of mn < 133 GeV is set at 95% confidence level. 
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1 Introduction 



The Standard Model (SM) of electroweak interactions 012] is tested with great precision by the 
experiments performed at the LEP and SLC e'^e" colhders running at centre-of-mass energies, 
^/s, close to the Z mass. From measurements of the total cross sections and forward-backward 
asymmetries in the reactions 

e^e~ hadrons(7) , e^e~ — > e^e~(7) , 

e"'"e~ /i"'"/i~(7) , e"'"e~ t^t~{^) , (1) 

the mass, total and partial widths of the Z and other electroweak parameters are obtained by 
L3 1311] and other experiments [SHHl- The (7) indicates the presence of radiative photons. 

The large luminosity collected in the years 1993 — 95 enables a significant improvement 
on our previous measurements of Z parameters. An integrated luminosity of 103 pb~^ was 
collected, corresponding to the selection of 2.5 ■ 10^ hadronic and 2.5 ■ 10^ leptonic events. 
Most of the data were collected at a centre-of-mass energy corresponding to the maximum 
annihilation cross section. 

In 1993 and 1995 scans, of the Z resonance were performed where runs at the Z pole 
alternated with runs at about 1.8 GeV on either side of the peak. Compared to previous 
measurements, our event samples on the wings of the Z resonance are increased by more than 
a factor of five. 

The LEP beam energies were precisely calibrated at the three energy points in 1993 — 95 
using the method of resonant depolarisation As a result, the contributions to the errors on 
the Z mass and total width from the uncertainty on the centre-of-mass energy are reduced by 
factors of about five and three, respectively, as compared to the data collected before. 

The installation of silicon strip detectors in front of the small angle electromagnetic calorime- 
ters allows a much more precise determination of the fiducial volume used for the luminosity 
measurement ^H]- This improvement, together with the reduced theoretical uncertainty on 
the small angle Bhabha cross section JJ^J^, allows more precise measurements of the cross 
sections, in particular that for e"^e~ — *■ hadrons(7). This results in a better determination of 
the invisible Z width, from which the number of light neutrino generations is deduced. 

In this article measurements of hadronic and leptonic cross sections and leptonic forward- 
backward asymmetries, obtained from the data collected between 1993 and 1995, are presented. 
These measurements are combined with our published results from the data collected in 1990 — 
92 j3]. The complete integrated luminosity collected by L3 at the Z resonance is 143 pb~^, 
consisting of about 3.5 ■ 10^ hadronic and 3.5 ■ 10^ leptonic events. The results on the properties 
of the Z boson and on other electroweak observables presented here are based on the final 
analyses of the complete data set collected at the Z resonance. 

This article is organised as follows: After a brief description of the L3 detector in Section |21 
we summarise in Section IHl features of the 1993 — 95 data analysis common to all final states 
investigated. Section |3] addresses issues related to the LEP centre-of-mass energy. The mea- 
surement of luminosity is described in Section |S| The event selection and the analysis of the 
reactions in (0) are discussed in Sections IHl to IHl and the results on the measurements of total 
cross sections and forward-backward asymmetries are presented in Section El A general de- 
scription of the fits performed to our data is given in Section ^2 Various fits for Z parameters 
are performed in Section El and the results of the fits in the framework of the SM are given in 
Section Uni We summarise and conclude in Section [T3| The Appendices 1X1 and iBl give details 
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on the treatment of the t-channel contributions in e^e — > e^e (7) and on technicahties of the 
fit procedures, respectively. 

2 The L3 Detector 

The L3 detector ^H] consists of a sihcon microvertex detector ^3], a central tracking chamber, 
a high resolution electromagnetic calorimeter composed of BGO crystals, a lead-scintillator ring 
calorimeter at low polar angles ^^1; a scintillation counter system, a uranium hadron calorime- 
ter with proportional wire chamber readout and an accurate muon spectrometer. Forward- 
backward muon chambers, completed for the 1995 data taking, extend the polar angle coverage 
of the muon system down to 24 degrees [TSj with respect to the beam line. All detectors are 
installed in a 12 m diameter magnet which provides a solenoidal field of 0.5 T in the central 
region and a toroidal field of 1.2 T in the forward-backward region. The luminosity is measured 
using BGO calorimeters preceded by silicon trackers [lOj situated on each side of the detector. 

In the L3 coordinate system the direction of the e~ beam defines the z direction. The xy^ 
or r0 plane, is the bending plane of the magnetic field, with the x direction pointing to the 
centre of the LEP ring. The coordinates and Q denote the azimuthal and polar angles. 

3 Data Analysis 

The data collected between 1993 and 1995 are split into nine samples according to the year 
and the centre-of-mass energy. Data samples at ^ ^ mz are referred to as peak, those at 
off-peak energies are referred to as peak— 2 and peak+2. The peak samples in 1993 and 1995 
are further split into data taken early in the year (pre-scan) and those peak runs interspersed 
with off-peak data taking (scan) which coincide with the precise LEP energy calibration (see 
Section Hj). Cross sections and leptonic forward-backward asymmetries are determined for each 
data sample. 

Acceptances, background contaminations and trigger efficiencies are studied for all nine data 
samples separately to take into account their possible dependence on the centre-of-mass energy 
and the time dependence of the detector status. Systematic errors are determined for the data 
samples individually. Average values for uncertainties are used if no dependence on the centre- 
of-mass energy or the data taking period is observed. Correlations of the systematic errors 
among the data sets are estimated and are taken into account in the analyses to determine 
electroweak parameters. 

Acceptances and background contaminations from e"^e~-interactions are determined by 
Monte Carlo simulations. The following event generator programs are used for the various 
signal and background processes: JETSET and HERWIG jIH] for e"'"e~ hadrons(7); 
KORALZ [19J for e+e" -> ;uV(7) and e+e"^^ r+r-(7); BHAGENE 120], BHWIDE [21] and 
BABAMC |22J for large angle e+e" ^ e+e" (7); BHLUMI [TJ for small angle e+e" ^ e+e" (7); 
GGG for e+e- ^ 77(7); DIAG36 |23j for e+e" ^ e+e'^+r; DIAG36, PHOJET ^3] and 
PYTHIA [r^ for e+e" — >■ e+e" hadrons. For the simulation of hadronic final states the frag- 
mentation parameters of JETSET and HERWIG are tuned to describe our data as discussed 
in Reference |26j . 

The generated events are passed through a complete detector simulation. The response of 
the L3 detector is modelled with the GEANT detector simulation program which includes 
the effects of energy loss, multiple scattering and showering in the detector materials. Hadronic 
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showers are simulated with the GHEISHA [2H] program. The performance of the detector, 
including inefficiencies and their time dependence as observed during data taking, is taken 
into account in the simulation. With this procedure, experimental systematic errors on cross 
sections and forward-backward asymmetries are minimized. 

4 LEP Energy Calibration 

The average centre-of-mass energy of the colliding particles at the L3 interaction point is calcu- 
lated using the results provided by the Working Group on LEP Energy [H]. Every 15 minutes 
the average centre-of-mass energy is determined from measured LEP machine parameters, ap- 
plying the energy model which is based on calibration by resonant depolarisation [20] • This 
model traces the time variation of the centre-of-mass energy of typically 1 MeV per hour. The 
average centre-of-mass energies are calculated for each data sample individually as luminosity 
weighted averages. Slightly different values are obtained for different reactions because of small 
differences in the usable luminosity. 

The errors on the centre-of-mass energies and their correlations for the 1994 data and for 
the two scans performed in 1993 and 1995 are given in form of a 7 x 7 covariance matrix in 
Table d The uncertainties on the centre-of-mass energy for the data samples not included in 
this matrix, i.e. the 1993 and 1995 pre-scans, are 18 MeV and 10 MeV, respectively. Details of 
the treatment of these errors in the fits can be found in Appendix iBl 

The energy distribution of the particles circulating in an e^e~-storage ring has a finite 
width due to synchrotron oscillations. An experimentally observed cross section is therefore a 
convolution of cross sections at energies which are distributed around the average value in a 
gaussian form. The spread of the centre-of-mass energy for the L3 interaction point as obtained 
from the observed longitudinal length of the particle bunches in LEP is listed in Table |2] [S]. 
The time variation of the average energy causes a similar, but smaller, effect which is included 
in these numbers. 

All cross sections and forward-backward asymmetries quoted below are corrected for the 
energy spread to the average value of the centre-of-mass energy. The relative corrections on 
the measured hadronic cross sections amount to +1.7 per mill (%c) at the Z pole and to — l.l%c 
and — 0.6%o at the peak— 2 and peak-|-2 energy, respectively. The absolute corrections on the 
forward-backward asymmetries are very small. The largest correction is —0.0002 for the muon 
and tau peak— 2 data sets. The error on the energy spread is propagated into the fits, resulting 
in very small contributions to the errors of the fitted parameters (see Appendix El- The largest 
effect is on the total width of the Z, contributing approximately 0.3 MeV to its error. 

During the operation of LEP, no evidence for an average longitudinal polarisation of the 
electrons or positrons has been observed. Stringent limits on residual polarisation during lumi- 
nosity runs are set such that the uncertainties on the determination of electroweak observables 
are negligible compared to their experimental errors jHU] . 

The determination of the LEP centre-of-mass energy in 1990 — 92 is described in Refer- 
ences jni]. From these results the LEP energy error matrix given in Table El is derived. 

5 Luminosity Measurement 

The integrated luminosity C is determined by measuring the number of small-angle Bhabha 
interactions e+e~ e+e~(7). For this purpose two cylindrical calorimeters consisting of arrays 
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of BGO crystals are located on either side of the interaction point. Both detectors are divided 
into two half-rings in the vertical plane to allow the opening of the detectors during filling of 
LEP. A silicon strip detector, consisting of two layers measuring the polar angle, 6', and one 
layer measuring the azimuthal angle, 0, is situated in front of each calorimeter to precisely 
define the fiducial volume. A detailed description of the luminosity monitor and the luminosity 
determination can be found in Reference [10]. 

The selection of small-angle Bhabha events is based on the energy depositions in adjacent 
crystals of the BGO calorimeters which are grouped to form clusters. The highest-energy 
cluster on each side is considered for the luminosity analysis. For about 98% of the cases a hit 
in the silicon detectors is matched with a cluster and its coordinate is used; otherwise the BGO 
coordinate is retained. 

The event selection criteria are: 

1. The energy of the most energetic cluster is required to exceed 0.8i?b and the energy on 
the opposite side must be greater than OAEy^, where Ei, is the beam energy. If the energy 
of the most energetic cluster is within ±5% of E-^, the minimum energy requirement on 
the opposite side is reduced to 0.2ii^b in order to recover events with energy lost in the 
gaps between crystals. The distributions of the energy of the most energetic cluster and 
the cluster on the opposite side as measured in the luminosity monitors are shown in 
Figure [U for the 1993 data. All selection cuts except the one under study are applied. 

2. The cluster on one side must be confined to a tight fiducial volume: 

• 32 mrad < ^ < 54 mrad; \cf> - 90°| > 11.25° and \4> - 270°| > 11.25°. 

The requirements on the azimuthal angle remove the regions where the half-rings of the 
detector meet. The cluster on the opposite side is required to be within a larger fiducial 
volume: 

• 27 mrad < n - 9 < 65 mrad; \(f) - 90°| > 3.75° and \4> - 270°| > 3.75°. 

This ensures that the event is fully contained in the detectors and edge effects in the 
reconstruction are avoided. 

3. The coplanarity angle A0 = (j){z < 0) — (f){z > 0) between the two clusters must satisfy 
|A0- 180°| < 10°. 

The distribution of the coplanarity angle is shown in Figure El Very good agreement with the 
Monte Carlo simulation is observed. 

Four samples of Bhabha events are defined by applying the tight fiducial volume cut to one 
of the 6'-measuring silicon layers. Taking the average of the luminosities obtained from these 
samples minimizes the effects of relative offsets between the interaction point and the detectors. 
The energy and coplanarity cuts reduce the background from random beam-gas coincidences. 
The remaining contamination is very small: (3.4 ± 2.2) ■ 10~^. This number is estimated using 
the sidebands of the coplanarity distribution, 10° < |A0 — 180° | < 30°, after requiring that 
neither of the two clusters have an energy within ±5% of E}^. 

The accepted cross section is determined from Monte Carlo e"''e~ — e"'"e~(7) samples gen- 
erated with the BHLUMI event generator at a fixed centre-of-mass energy of ^/s = 91.25 GeV. 
The dependence on the centre-of-mass energy, as well as the contributions of Z-exchange and 
7Z interference, are calculated with the BHLUMI program. At y/s = 91.25 GeV the accepted 
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cross section is determined to be 69.62 nb. The statistical error on the Monte Carlo sample con- 
tributes 0.35%c to the uncertainty of the luminosity measurement. The theoretical uncertainty 
on the Bhabha cross section in our fiducial volume is estimated to be 0.61%o |12j . 

The experimental errors of the luminosity measurement are small. Important sources of 
systematic errors are: geometrical uncertainties due to the internal alignment of the silicon 
detectors (0.15%e to 0.27%e), temperature expansion effects (0.14%c) and the knowledge on the 
longitudinal position of the silicon detectors (0.16%e to 0.60%c). The precision depends on the 
accuracy of the detector surveys and on the stability of the detector and wafer positions during 
the different years. 

The polar angle distribution of Bhabha scattering events used for the luminosity measure- 
ment is shown in Figure El The structure seen in the central part of the +z side is due to the 
flare in the beam pipe on this side. The imperfect description in the Monte Carlo does not pose 
any problem as it is far away from the edges of the fiducial volume. 

The overall agreement between the data and Monte Carlo distributions of the selection 
quantities is good. Small discrepancies in the energy distributions at high energies are due to 
contamination of Bhabha events with beam-gas interactions and, at low energies, due to an 
imperfect description of the cracks between crystals. The selection uncertainty is estimated by 
varying the selection criteria over reasonable ranges and summing in quadrature the resulting 
contributions. This procedure yields errors between 0.42%o and 0.48%c for different years. The 
luminosities determined from the four samples described above agree within these errors. The 
trigger inefficiency is measured using a sample of events triggered by only requiring an energy 
deposit exceeding 30 GeV on one side. It is found to be negligible. 

The various sources of uncertainties are summarized in Table El Combining them in quadra- 
ture yields total experimental errors on the luminosity of 0.86%c, 0.64%c and 0.68%o in 1993, 1994 
and 1995. Correlations of the total experimental systematic errors between different years are 
studied and the correlation matrix is given in Table The error from the theory is fully 
correlated. 

Because of the 1/ s dependence of the small angle Bhabha cross section, the uncertainty on 
the centre-of-mass energies causes a small additional uncertainty on the luminosity measure- 
ment. For instance, this amounts to 0.1%c for the high statistics data sample of 1994. This 
effect is included in the fits performed in Section [T21 and see Appendix |Bl 

The statistical error on the luminosity measurement from the number of observed small 
angle Bhabha events is also included in those fits. Table (HI lists the number of observed Bhabha 
events for the nine data samples and the corresponding errors on cross section measurements. 
Combining all data sets taken in 1993 — 95 at ^/s ~ rrtz the statistical error on the luminosity 
contributes 0.45%c to the uncertainty on the pole cross section measurements. 

Higher order corrections from photon radiation to the small angle Bhabha cross section 
are studied with the photon spectrum of luminosity events. For this analysis events with two 
distinct energy clusters exceeding O.li^b in one of the calorimeters are selected. The photon 
is identified as the lower energy cluster. The fraction of radiative events with > O.li^b in 
the total low-angle Bhabha sample is 2% and the measured cross section, normalised to the 
expectation, is found to be 0.993 ±0.16. The observed spectrum from 1993 is shown in FigureEl 
and good agreement is found with the Monte Carlo expectation. 
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6 e+e — > hadrons(7) 

Event Selection 

Hadronic Z decays are identified by their large energy deposition and high multiplicity in the 
electromagnetic and hadron calorimeters. The selection criteria are similar to those applied in 
our previous analysis 0: 

1. The total energy observed in the detector, E^is, normalised to the centre-of-mass energy 
must satisfy 0.5 < E^is/^/s < 2.0; 

2. The energy imbalance along the beam direction, E\\, must satisfy |/i5^vis < 0.6; 

3. The transverse energy imbalance, E±, must satisfy E±/E^is < 0.6; 

4. The number of clusters, A^ci, formed from energy depositions in the calorimeters is required 
to be: 

a) Nci > 13 for | cos6t\ < 0.74 (barrel region), 

b) A'ci > 17 for I cos6't| > 0.74 (end-cap region), 
where 6t is the polar angle of the event thrust axis. 

Detailed analyses of the large data samples collected have been used to improve the Monte 
Carlo simulation of the detector response. FiguresElto|nishow the distributions of the quantities 
used to select hadronic Z decays and the comparisons to the Monte Carlo predictions. In these 
plots all selection cuts are applied, except the one under study. Good agreement is observed 
between our data and the Monte Carlo simulations. 

Total Cross Section 

The acceptance for e+e^ —>■ hadrons(7) events is determined from large samples of Monte Carlo 
events generated with the JETSET program. Applying the selection cuts, between 99.30% and 
99.42% of the events are accepted depending on the year of the data taking and on differences 
in initial-state photon radiation at the various centre-of-mass energies. Monte Carlo events are 
generated with > O.ly^ where y/s' is the effective centre-of-mass energy after initial state 
photon radiation. The acceptance for events in the data with \/s' < O.ly^ is estimated to be 
negligible. They are not considered as part of the signal and hence not corrected for. 

The interference between initial and final state photon radiation is not accounted for in the 
event generator. This effect modifies the angular distribution of the events in particular at 
very low polar angles where the detector inefficiencies are largest. However, the error from the 
imperfect simulation on the measured cross section, which includes initial-final state interference 
as part of the signal, is estimated to be very small (<C 0.1 pb) in the centre-of-mass energy 
range considered here. Quark pairs originating from pair production from initial state radiation 
are considered as part of the signal if their invariant mass exceeds 50% of ^/s. 

To estimate the uncertainty on the acceptance on the modelling of the quark fragmentation, 
the determination of the acceptance is repeated using the HERWIG program. The detector 
simulations of both Monte Carlo programs are tuned in the same way to describe as closely as 
possible our data, e.g. in terms of energy resolution and cluster multiplicity. The remaining 
difference in acceptance is 0.42%c and we assign half of it as an estimate of the uncertainty 
on the acceptance of e"'"e~ —>■ hadrons(7) events due to the modelling of quark fragmentation. 
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Differences of the implementation of QED effects in both programs are studied and found to 
have neghgible impact on the acceptance. 

Hadronic Z decays are triggered by the energy, central track, muon or scintillation counter 
multiplicity triggers. The combined trigger efficiency is obtained from the fraction of events 
with one of these triggers missing as a function of the polar angle of the event thrust axis. 
This takes into account most of the correlations among triggers. A sizeable inefficiency is only 
observed for events in the very forward region of the detector, where hadrons can escape through 
the beam pipe. Trigger efficiencies, including all steps of the trigger system, between 99.829% 
and 99.918% are obtained for the various data sets. Trigger inefficiencies determined for data 
sets taken in the same year are statistically compatible. Combining those data sets results in 
statistical errors of at most 0.12%c which is assigned as systematic error to all data sets. 

The background from other Z decays is found to be small: 2.9%o essentially only from 
e"'"e~ T^T~ {'-/). The uncertainty on this number is negligible compared to the total systematic 
error. 

The determination of the non-resonant background, mainly e"'"e~ — >■ e"'"e~ hadrons, is based 
on the measured distribution of the visible energy shown in Figure The Monte Carlo program 
PHOJET is used to simulate two-photon collision processes. The absolute cross section is 
derived by scaling the Monte Carlo to obtain the best agreement with our data in the low end 
of the ii^vis spectrum: 0.32 < E^^s/ y/s < 0.44. Consistently for all data sets, scale factors of 1.1 
are necessary. In the signal region contaminations from e^e~ — >■ e+e~ hadrons between 11.6 pb 
and 13.0 pb are obtained for the different data sets. No dependence on ^/s is observed. This is 
in agreement with results of a similar calculation performed with the DIAG36 program. 

Beam related background (beam-gas and beam-wall interactions) is small. To the extent 
that the E^is spectrum is similar to that of e+e~ —>■ e"'"e~ hadrons, it is accounted for by 
determining the absolute normalisation from the data. 

As a check, the non-resonant background is estimated by extrapolating an exponential 
dependence of the i?vis spectrum from the low energy part into the signal region. This method 
yields consistent results. Based on these studies we assign an error on the measured hadron 
cross section of 3 pb due to the understanding of the non-resonant background. This error 
assignment is supported by our measurements of the hadronic cross section at high energies 
(130 GeV < y/s < 172 GeV) where the relative contribution of two-photon processes is much 
larger [^1221 • The extrapolation of these studies back to the Z peak yields a similar result for 
the uncertainty. 

The contribution of random uranium noise and electronic noise in the detector faking a 
signal event is determined from a subsample of the event candidates. This subsample is obtained 
requiring that most of the observed energy stems either from the electromagnetic or the hadron 
calorimeter and that there be little matching between individual energy deposits and tracks. 
The Ej_/E^is distribution of this subsample shows an e"'"e~ — >• hadrons(7) signal over a fiat 
background (see Figure HO for the 1994 data). This background is consistent with a constant 
noise rate, from which a background correction of 7.4 pb is derived. An uncertainty of 1 pb 
on the hadron cross section is assigned to all data sets from this correction. The absolute 
normalisation of the e"'"e~ hadrons(7) signal in Figure ^1 is not expected to be perfectly 
reproduced by the Monte Carlo simulation. However, this does not pose a serious problem as 
the noise rate is determined from the tail of the spectrum. 

The systematic error from event selection on the measured cross sections is estimated by 
varying the selection cuts. All cross section results are stable within ±0.3%c. The systematic 
errors to the cross section measurements e"'"e~ —>■ hadrons(7) are summarised in Table [7| Un- 
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certainties which scale with the cross section and absolute uncertainties are separated because 
they translate in a different way into errors on Z parameters, in particular on the total width. 
The scale error is further split into a part uncorrelated among the data samples, in this case 
consisting of the contribution of Monte Carlo statistics, and the rest which is taken to be fully 
correlated and amounts to 0.39%o. 

The results of the e"'"e~ — > hadrons(7) cross section measurements are discussed in Sec- 
tion EB 

7 e+e~ /x+/x~(7) 

Event Selection 

The selection of e"'"e~ in the 1993 and 1994 data is similar to the selection applied in 

previous years described in Reference Two muons in the polar angular region | cos^| < 0.8 
are required. Most of the muons, 88%, are identified by a reconstructed track in the muon 
spectrometer. Muons are also identified by their minimum ionising particle (MIP) signature 
in the inner sub-detectors, if less than two muon chamber layers are hit. A muon candidate is 
denoted as a MIP, if at least one of the following conditions is fulfilled: 

1. A track in the central tracking chamber must point within 5° in azimuth to a cluster in 
the electromagnetic calorimeter with an energy less than 2 GeV. 

2. On a road from the vertex through the barrel hadron calorimeter, at least five out of a 
maximum of 32 cells must be hit, with an average energy of less than 0.4 GeV per cell. 

3. A track in the central chamber or a low energy electromagnetic cluster must point within 
10° in azimuth to a muon chamber hit. 

In addition, both the electromagnetic and the hadronic energy in a cone of 12° half-opening 
angle around the MIP candidate, corrected for the energy loss of the particle, must be less than 
5 GeV. 

Events of the reaction e"^e" fi^fi~{'-f) are selected by the following criteria: 

1. The event must have a low multiplicity in the calorimeters Nd < 15. 

2. If at least one muon is reconstructed in the muon chambers, the maximum muon momen- 
tum must satisfy Pmax > 0.6 E]^. If both muons are identified by their MIP signature there 
must be two tracks in the central tracking chamber with at least one with a transverse 
momentum larger than 3 GeV. 

3. The acollinearity angle must be less than 90°, 40° or 5° if two, one or no muons are 
reconstructed in the muon chambers. 

4. The event must be consistent with an origin of an e"'"e~-interaction requiring at least one 
time measurement of a scintillation counter, associated to a muon candidate, to coincide 
within ±3 ns with the beam crossing. Also, there must be a track in the central tracking 
chamber with a distance of closest approach to the beam axis of less than 5 mm. 

As an example. Figure ITD shows the distribution of the maximum measured muon momen- 
tum for candidates in the 1993 — 94 data compared to the expectation for signal and background 
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processes. The acoUinearity angle distribution of the selected muon pairs is shown in Figure IT2l 
The experimental angular resolution and radiation effects are well reproduced by the Monte 
Carlo simulation. 

The analysis of the 1995 data in addition uses the newly installed forward-backward muon 
chambers. The fiducial volume is extended to |cos^| < 0.9. Each event must have at least 
one track in the central tracking chamber with a distance of closest approach in the transverse 
plane of less than 1 mm and a scintillation counter time coinciding within ±5 ns with the beam 
crossing. The rejection of cosmic ray muons in the 1995 data is illustrated in Figure El 

For events with muons reconstructed in the muon chambers the maximum muon momentum 
must be larger than | i?b. Every muon without a reconstructed track in the muon chambers 
must have a transverse momentum larger than 3 GeV as measured in the central tracking 
chamber. The polar angle distribution of muon pairs collected in 1995 is shown in Figure El 

Total Cross Section 

The acceptance for the process e"'"e~ fi~^^~{j) in the fiducial volume | cos6'| < 0.8 (0.9 for 
1995 data) and for C, < 90° is determined with events generated with the KORALZ program. 
We obtain acceptances between 92.25% and 93.04%, mainly depending on the centre-of-mass 
energy. The systematic error on the cross section from imperfect description of detector inef- 
ficiencies is estimated to be 2.7%c (3.2%o for the 1995 data). This number is calculated from a 
comparison with results obtained by removing events at the detector edges from the analysis 
and using different descriptions of time dependent detector inefficiencies. Smaller contribu- 
tions to the systematic error arise from the statistical precision of the Monte Carlo simulations 
performed for the different data samples. 

Muon pairs are mainly triggered by the muon and the central track trigger. The trigger 
efficiencies are studied as a function of the azimuthal angle as inefficiencies are expected close 
to chamber boundaries. For the 1995 data also the polar angular dependence of the trigger 
efficiency is determined to account for effects in the forward region. Events with both muons 
reconstructed in the muon chambers are triggered with full efficiency. The efficiency of the 
central track trigger is independently determined using Bhabha events. The overall trigger 
efficiency varies between 99.62% and 99.90% for the different years of data taking. Systematic 
errors on the measured cross sections of less than l%c are estimated from comparing a simulation 
of the central track trigger efficiency and its measurement with Bhabha events. 

A background of (1.35 ± 0.03)% remains in the sample arising from e^e~ r+r~(7) 
events with both tau leptons decaying into muons. The error reflects Monte Carlo statistics 
and the uncertainty of the branching ratio yU~z/^z/^ [34 . Other backgrounds from Z 

decays are smaller than 0.1%o. The contamination from the non- resonant two-photon process 
e"'"e~ —>■ e'^e~ n'^fi~ is 0.11 pb, i.e. between 0.1%o and 0.3%c of the signal cross section, as 
determined using the DIAG36 Monte Carlo program. 

The residual contamination from cosmic ray muons in the event sample is determined from 
the sideband in the distribution of distance of closest approach to the beam axis after all other 
selection cuts are applied (Figure EJ- Cosmic ray muons enter into the event sample at a rate 
of (9.7 ± 0.8) ■ 10~^ per minute of data taking which translates to background contaminations 
between 1.9%e and 6.8%c for the different data sets depending on their average instantaneous 
luminosity and the signal cross section. The statistical precision of the determination of the 
cosmic contamination causes a systematic error of 0.3 pb on the total muon pair cross section. 

By varying the selection cuts we determine systematic errors on the total cross section 
between 1.3%o and 2.2%o. The systematic errors on the cross section measurements e"'"e~ — *■ 
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/i'^/i~(7) are summarised in Table |H1 

Resonant four-fermion final states with a high-mass muon pair and a low-mass fermion 
pair are accepted. These events are considered as part of the signal if the invariant mass of 
the muon pair exceeds 0.5^/s. This inclusive selection minimizes errors due to higher order 
radiative corrections. Especially no cut is applied on additional tracks from low-mass fermion 
pairs in the final state . 



Forward-Backward Asymmetry 

The forward-backward asymmetry, ApB, is defined as: 

A™ = ^, (2) 

CTp + CTb 

where ap is the cross section for events with the fermion scattered into the hemisphere which is 
forward with respect to the e~ beam direction. The cross section in the backward hemisphere 
is denoted by o"b. Events with hard photon bremsstrahlung are removed from the sample 
by requiring that the acollinearity angle of the event be less than 15°. The differential cross 
section in the angular region | cos6'| < 0.9 can then be approximated by the lowest order angular 
dependence to sufficient precision: 

da 3 



oc 



dcos^ 8 



(l + cos^^) +Afbcos^, (3) 



with 6 being the polar angle of the final state fermion with respect to the e~ beam direction. 

For each data set the forward-backward asymmetry is determined from a maximum likeli- 
hood fit to our data where the likelihood function is defined as the product over the selected 
events labelled i of the differential cross section evaluated at their respective scattering angle 

e^■. 

^ = n (l + cos' Oi) + (1 - 2k,) Afb cos . (4) 

The probability of charge confusion for a specific event, Kj, is included in the fit. Only events 
with opposite charge assignment to the two muons are used for this measurement. The bias 
on the asymmetry measurement introduced by the use of the lowest order angular dependence 
(Equation 01) does not exceed 0.0003. 

This method does not require an exact knowledge of the acceptance as a function of the 
polar angle provided that the acceptance is independent of the muon charge. Events without a 
reconstructed muon in the muon chambers are included with the charge assignment obtained 
from the central tracking chamber in a similar way as for e"'"e~ final states This largely 
reduces effects of charge dependent acceptance in the muon chambers. The remaining asymme- 
try is estimated by artificially symmetrising the detector. For each known, inefficient detector 
element, the element opposite with respect to the centre of the detector is removed from the 
data reconstruction. The event selection is applied again and, for the large 1994 data set, the 
measured forward-backward asymmetry changes by 0.0011 ± 0.0006. Half of this difference, 
0.0006, is assigned to all data sets as a systematic error on ApB from a possible detector asym- 
metry. In 1995 the forward-backward muon chambers did not contribute significantly to the 
detector asymmetry. 

The values of Ki are obtained from the fraction of events with identical charges assigned to 
both muons. Besides its dependence on the transverse momentum, the charge measurement 
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strongly depends on the number of muon chamber layers used in the reconstruction. The 
charge confusion is determined for each event class individually. The average charge confusion 
probability, almost entirely caused by muons only measured in the central tracking chamber, 
is (3.2 ± 0.3)%o, (0.8 ± 0.1)%o and (1.0 ± 0.3)%c for the years 1993, 1994 and 1995, respectively, 
where the errors are statistical. The improvement in the charge determination for 1994 and 
1995 reflects the use of the silicon microvertex detector. 

The correction for charge confusion is proportional to the forward-backward asymmetry and 
it is less than 0.001 for all data sets. To estimate a possible bias from a preferred orientation 
of events with the two muons measured to have the same charge we determine the forward- 
backward asymmetry of these events using the track with a measured momentum closer to the 
beam energy. The asymmetry of this subsample is statistically consistent with the standard 
measurement. Including these like-sign events in the 1994 sample would change the measured 
asymmetry by 0.0008. Half of this number is taken as an estimate of a possible bias of the 
asymmetry measurement from charge confusion in the 1993 — 94 data. The same procedure is 
applied to the 1995 data and the statistical precision limits a possible bias to 0.0010. 

Differences of the momentum reconstruction in forward and backward events would cause 
a bias of the asymmetry measurement because of the requirement on the maximum measured 
muon momentum. We determine the loss of efficiency due to this cut separately for forward and 
backward events by selecting muon pairs without cuts on the reconstructed momentum. No 
significant difference is observed and the statistical error of this comparison limits the possible 
effect on the forward-backward asymmetry to be less than 0.0004 and 0.0009 for the 1993 — 94 
and 1995 data, respectively. 

Other possible biases from the selection cuts on the measurement of the forward-backward 
asymmetry are negligible. This is verified by a Monte Carlo study which shows that events not 
selected for the asymmetry measurement, but inside the fiducial volume and with ^ < 15°, do 
not have a different Ap^ value. 

The background from e"'"e~ r^r^(7) events is found to have the same asymmetry as the 
signal and thus neither necessitates a correction nor causes a systematic uncertainty. The effect 
of the contribution from the two-photon process e"'"e~ — >• e"'"e~ /i"*"/!", further reduced by the 
tighter acollinearity cut on the measured muon pair asymmetry, can be neglected. The forward- 
backward asymmetry of the cosmic ray muon background is measured to be —0.02 ±0.13 using 
the events in the sideband of the distribution of closest approach to the interaction point. 
Weighted by the relative contribution to the data set this leads to corrections of —0.0007 and 
-1-0.0003 to the peak— 2 and peak-|-2 asymmetries, respectively. On the peak this correction is 
negligible. The statistical uncertainty of the measurement of the cosmic ray asymmetry causes 
a systematic error of 0.0001 on the peak and between 0.0003 and 0.0005 for the peak— 2 and 
peak-|-2 data sets. 

The systematic uncertainties on the measurement of the muon forward-backward asymmetry 
are summarised in Table IHl In 1993 — 94 the total systematic error amounts to 0.0008 at the 
peak points and to 0.0009 at the off-peak points due to the larger contamination of cosmic ray 
muons. For the 1995 data the determination of systematic errors is limited by the number of 
events taken with the new detector configuration and the total error is estimated to be 0.0015. 

In Figure IT5l the differential cross sections dcr/dcos^ measured from the 1993 — 95 data 
sets are shown for three different centre-of-mass energies. The data are corrected for detector 
acceptance and charge confusion. Data sets with a centre-of-mass energy close to mz, as well 
as the data at peak— 2 and the data at peak -(-2, are combined. The data are compared to the 
differential cross section shape given in Equation |S| 
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The results of the total cross section and forward-backward asymmetry measurements in 
e"'"e~ — ^ /i+/i~(7) are presented in Section ITHl 

8 e+e~ r+T~(7) 

Event Selection 

The selection of e"'"e~ r+r~(7) events aims to select all hadronic and leptonic decay modes 
of the tau. Z decays into tau leptons are distinguished from other Z decays by the lower visible 
energy due to the presence of neutrinos and the lower particle multiplicity as compared to 
hadronic Z decays. Compared to our previous analysis ^ the selection of e"'"e~ t~^t~{'j) 
events is extended to a larger polar angular range, | cos6'i| < 0.92, where 9t is defined by the 
thrust axis of the event. 

Event candidates are required to have a jet, constructed from calorimetric energy de- 
posits fHU] and muon tracks, with an energy of at least 8 GeV. Energy deposits in the hemisphere 
opposite to the direction of this most energetic jet are combined to form a second jet. The two 
jets must have an acollinearity angle < 10°. There is no energy requirement on the second 
jet. 

High multiplicity hadronic Z decays are rejected by allowing at most three tracks matched 
to any of the two jets. In each of the two event hemispheres there should be no track with an 
angle larger than 18° with respect to the jet axis. Resonant four-fermion final states with a high 
mass tau pair and a low mass fermion pair are mostly kept in the sample. The multiplicity cut 
affects only tau decays into three charged particles with the soft fermion close in space leading 
to corrections of less than l%c. 

If the energy in the electromagnetic calorimeter of the first jet exceeds 85%, or the energy of 
the second jet exceeds 80%, of the beam energy with a shape compatible with an electromagnetic 
shower the event is classified as e"'"e~ — >• e"'"e~(7) background and hence rejected. 

Background from e"'"e~ — > fi'^fi~ (7) is removed by requiring that there be no isolated muon 
with a momentum larger than 80% of the beam energy and that the sum of all muon momenta 
does not exceed 1.5E]^. Events are rejected if they are consistent with the signature of two 
MIPs. 

To suppress background from cosmic ray events the time of scintillation counter hits asso- 
ciated to muon candidates must be within ±5 ns of the beam crossing. In addition, the track 
in the muon chambers must be consistent with originating from the interaction point. 

In Figures El to El the energy in the most energetic jet, the number of tracks associated 
to both jets, the acollinearity between the two jets and the distribution of | cos6'(| are shown 
for the 1994 data. Data and Monte Carlo expectations are compared after all cuts are applied, 
except the one under study. Good agreement between data and Monte Carlo is observed. Small 
discrepancies seen in Figure IT7I are due to the imperfect description of the track reconstruction 
efficiency in the central chamber. Their impact on the total cross section measurement is small 
and is included in the systematic error given below. 

Tighter selection cuts must be applied in the region between barrel and end-cap part of 
the BGO calorimeter and in the end-cap itself, reducing the selection efficiency (see Figure El) • 
This is due to the increasing background from Bhabha scattering. Most importantly the shower 
shape in the hadron calorimeter is also used to identify candidate electrons and the cuts on the 
energy of the first and second jet in the electromagnetic end-cap calorimeter are tightened to 
75% of the beam energy. 
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Total Cross Section 

Between 70.21% and 70.91% of the signal events are accepted inside the fiducial volume de- 
fined by I cos 9t\ < 0.92. The acceptance for e"'"e~ r+r~(7) events depends on the tau decay 
products. The experimental knowledge of tau branching fractions jSl] translates to an uncer- 
tainty on the average acceptance of e"'"e~ r+r^(7) events which contributes with 2%o to the 
systematic error on the cross section measurement. From the data the efficiency of the trigger 
system for selected e"'"e~ — >• t~^t~{'j) events is determined to be (99.71 ± 0.02)%. 

The largest remaining background consists of Bhabha events, 1.3% to 3.7%, depending 
on the centre-of-mass energy, entering into the sample predominantly at low polar angles. 
Background from Z decays into hadrons is determined to be between 1.3%c and 2.7%c, depending 
on the data taking period, and 7.5%c from Z decays into muons. The statistical precision of 
the background determination by Monte Carlo simulations causes systematic errors between 
1.0%o and 3.3%e. Contaminations from non-resonant background are small: l%c to 2%o from 
two-photon collisions and 2%c to 3%c from cosmic ray muons, depending on the centre-of-mass 
energy. The systematic error from the subtraction of non-resonant background is estimated to 
be 1.2 pb. 

From variations of the above selection cuts contributions to the systematic error on the total 
cross section between 5.3%c and 8.0%c are estimated for different years, largely independent of 
the centre-of-mass energy. The main contribution arises from the definition of the fiducial 
volume by | cos 6t\ < 0.92, see Figure IT^ The systematic errors on the e^e^ r+r^(7) cross 
section measurements are summarised in Table ITOl 

Forward-Backward Asymmetry 

The forward-backward asymmetry of e^e^ r+r^(7) events is determined in the same way 
as described for muon pairs (Equation |^. The charge of a tau is derived from the sum of the 
charges of its decay products as measured in the central tracking and the muon chambers. The 
event sample selected for the cross section measurement is used requiring opposite and unit 
charge for the two tau jets. 

The average probability for a mis-assignment of both charges as determined from the ratio 
of like and unlike sign events is (7.4±0.4)%o in 1993. The use of the silicon microvertex detector 
reduced this mis-assignment to (2.5 ± 0.1)%c and (1.3 ± 0.1)%c in 1994 and 1995. Because the 
charge confusion probability is approximately independent of the polar angle this average value 
is used in the fit for A-pB- The systematic error on the forward-backward asymmetry from 
the uncertainty in the determination and the treatment of the charge confusion probability is 
estimated to be less than 0.0001 for all data sets. 

The effect of a possible detector asymmetry, in particular at the edges of the fiducial volume, 
is estimated from variation of the cos 6t cut. The statistical accuracy of this test limits this 
uncertainty to 0.003 which is taken as a systematic error. The measured asymmetries are 
corrected for background contributions. The uncertainty on the background contamination, in 
particular from e^e~ — >■ e"'"e~(7), translates into an error of 0.001 on the tau pair asymmetry. 

Large Monte Carlo samples are used to study a possible bias on the measured asymme- 
try from the fit method and from the selection cuts. In particular, energy and momentum 
requirements might preferentially select certain helicity configurations leading to a bias in the 
determination of Afb- The Monte Carlo simulation does not show evidence for such a bias and 
its statistical precision, 0.0004, is taken as the systematic error. 

During the 1995 data taking, large shifts of the longitudinal position of the e+e~-interaction 
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point were observed caused by the reconfiguration of tlie LEP radio frequency system [Oj. How- 
ever, tliey are found to have no sizeable effect on the measurement of the forward-backward 
asymmetry. The total systematic error assigned to the forward-backward asymmetry measure- 
ment of tau pairs is 0.0032 fTablelTT|). It is fully correlated between the data sets. 

The measured differential cross sections, combining the data into three centre-of-mass energy 
points, are shown in Figure l20l The lines show the results of fits to the data using the functional 
form of Equation IHl 

Section El presents the measurements of the total cross section and the forward-backward 
asymmetry in e"'"e~ T'^T~{'y). 

9 e+e~ — > e+e~(7) 

Event Selection 

The analysis of the reaction e"'"e~ e+e~ (7) is restricted to the polar angular range 44° < 6 < 
136° to increase the relative contribution of Z exchange to the measured cross section. The sig- 
nature of e"'"e~ final states is the low multiplicity high energy deposition in the electromagnetic 
calorimeter with associated tracks in the central tracking chamber. 

Most of the events are selected by requiring at least two clusters in the fiducial volume of 
the electromagnetic calorimeter, one with an energy greater than 0.9 i?b and the other with 
more than 2 GeV. The polar angles are determined form the centre-of-gravity of the clusters 
in the calorimeter and the interaction point. Figure 1^ shows the distribution of the highest 
energy cluster, Ei, normalised to the beam energy for events which pass all cuts except the 
requirement on the most energetic cluster. 

Electrons are discriminated from photons by requiring five out of 62 anodes of the central 
tracking chamber with a hit matching in azimuthal angle within ±3° with the cluster in the 
calorimeter. Two electron candidates are required inside the fiducial volume and with an 
acoUinearity angle ^ < 25°. Figure |221 shows the distribution of the acoUinearity angle. All 
other cuts except the one under study are applied. 

The event selection depends on the exact knowledge of imperfections of the electromag- 
netic calorimeter. The impact of the discrepancies seen in Figure |^ around the cut value is 
significantly reduced by accepting also events without a second cluster in the fiducial volume 
of the electromagnetic calorimeter. In this case a cluster in the hadron calorimeter is required 
consistent with an electromagnetic shower shape and at least 7.5 GeV opposite to the leading 
BGO cluster. This recovers events, up to 4%c of the total sample, with electrons leaking through 
the BGO support structure. Events failing the requirement on the most energetic cluster in the 
electromagnetic calorimeter are accepted if the sum of the energies of the four highest energy 
clusters anywhere in the electromagnetic calorimeter is larger than 70% of the centre-of-mass 
energy. In addition this partially recovers radiative events. 

For all event candidates the total number of energy deposits, A^^d, must be less than 15 (12 
for 1995 data). 

Total Cross Section 

The selection efficiency is determined using Monte Carlo events generated with the program 
BHAGENE, which generates up to three photons. Efficiencies between 97.37% and 98.53% are 
obtained for the different samples, where the differences originate from time dependent detector 
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inefficiencies. The use of high multiplicity hadron events allows to monitor the status of each 
individual BGO crystal in short time intervals. Inefficient crystals, typically 100 out of 8000 
in the barrel part, are identified and taken into account in the Monte Carlo simulation. This 

method, together with the redundancy of the selection cuts, reduces the systematic error on 
the selection efficiency. Limited Monte Carlo statistics causes systematic errors between 0.4%c 
and 1.0%c. 

The calculation of the selection efficiency is checked using events generated with the pro- 
grams BABAMC and BHWIDE. The efficiencies calculated with the different event generators 
agree within ±l%o which is taken as an estimate of the systematic error. 

The efficiency of the electron and photon discrimination in the central tracking chamber 
is determined using a subsample of data events selected by a tight acoUinearity cut < 1°) 
and requiring two high energy clusters in the electromagnetic calorimeter (E > 30 GeV) . Here 
the contamination of c+c" e"'"e~7 events with one electron and the photon inside, and the 
other electron outside the fiducial volume is expected to be very small. In this sample, events 
with only one identified electron originate from mis-identified Bhabha events or from photon 
conversion of e+e~ — > 77 events. The contamination of the latter in this sample is 0.4%c to 0.9%c 
as calculated from Monte Carlo. After correction for this contamination, the probability that 
one of the electrons in e"'"e~ final states fails the electron-photon discrimination is measured to 
be (0.7 ± 0.3)%c and (1.0 ± 0.1)%c for the 1993 and 1994 data, respectively. We correct for this 
effect. 

The method to determine this probabihty from the data is checked on fully simulated 
e'^e" — > e"^e~(7) Monte Carlo events. Firstly by not applying the electron- photon discrimina- 
tion, the contamination of events in the data used for the cross section measurement with one 
photon and only one electron in the fiducial volume is determined to be (2.6 ± 0.5)%o. This is 
in reasonable agreement with the Monte Carlo prediction of 1.4%c. Then we apply the above 
method to determine the probability that an electron fails the electron-photon discrimination 
on the fully simulated events and compare it to the value obtained using the generator informa- 
tion. The result is consistent within 0.6%o which is assigned as a systematic error to the total 
cross section due to the simulation and determination of the electron-photon discrimination. 

In 1995 the quality criteria on the status of the central tracking chamber are relaxed to 
increase the data sample at the expense of a smaller efficiency on the electron identification 
and a larger systematic error. Between 1.9%c and 2.8%c of the electrons fail the electron- photon 
discrimination cuts as determined from Monte Carlo simulation. We correct for this effect and 
a systematic error of 1.5%o is assigned to the total cross section measurement. 

Large angle Bhabha scattering events are triggered by the energy and the central track 
triggers. The overall trigger inefficiency is found to be < 0.1%c and has a negligible effect on 
the cross section measurement. 

In the 1993 and 1994 data the longitudinal position of the e"'"e~ interaction point is stable 
within ±2 mm. The corresponding uncertainty on the definition of the fiducial volume translates 
to a systematic error of 0.5%c on the cross section measurement. Imperfections of the description 
of the BGO geometry and the shower shape of electrons lead to a possible difference of the 
definition of the polar angle between data and Monte Carlo simulation. This difference is found 
to be less than 0.1°, translating to a systematic error of 0.5%o on the cross section measurement. 

The large movements of the interaction point in 1995 are determined from our e"'"e~ — >■ 
hadrons(7) data and the positions are used to calculate the scattering angle in e"'"e~ — * e"'"e~(7) 
events. The remaining systematic uncertainty on the definition of the fiducial volume, including 
the description of the BGO geometry, is estimated from a variation of the cut on the polar angle 
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to be 1.5%o. 

The selected sample contains about 1% background from the process e^e^ r^r^ {'-/), only 
slightly depending on the centre-of-mass energy. Contaminations from hadronic Z decays and 
the process e"'"e~ — >■ e"'"e~ e^e~ are below l%o and the remaining background from e'^e" — >■ 77 
is negligible. The error on the total cross section from background subtraction is 0.4%c to 1.0%o 
originating from limited Monte Carlo statistics. 

The systematic uncertainty of the event selection, estimated from variations of the selection 
cuts around their nominal values, varies between 0.8%o and 2.7%c for the various data sets. The 
systematic uncertainties contributing to the measurement of the cross section e"'"e~ — >■ e"'"e~(7) 
are summarised in Table IT^ 



Forward-Backward Asymmetry 

The data sample for the forward-backward asymmetry measurement is obtained from the sam- 
ple used for the measurement of the total cross section requiring in addition that each of the 
two electron candidates match with a track within 25 mrad in azimuth. 

The charge determination of the electrons is described in detail in reference The charge 
confusion is measured with the e"'"e~ —>■ (7) data sample which has an independent charge 
measurement from the muon spectrometer. We obtain for the probability of a wrong event 
orientation values between 0.5% and 4.6%. Lower values are due to the exploitation of the 
silicon microvertex detector in 1994 and 1995. We determine the asymmetry of a subsample 
with much lower charge confusion by excluding events with tracks close to the cathode and 
anode planes of the central tracking chamber. Comparing these results to those obtained from 
the full sample we derive a systematic error on Ap^ of 0.002 from the uncertainty of the charge 
determination. 

In the event sample used for the asymmetry measurement the main background from 
e"'"e~ r+r^(7) is reduced to about 4%c because the tight requirement on the matching 
between tracks and clusters in the electromagnetic calorimeter removes r~ — ^ P~i^t decays 
present in the cross section sample. It induces a correction of 0.002 on the asymmetry for the 
peak— 2 and of less than 0.0005 for the other data sets. The effect is largest at peak— 2 because 
of the difference of the e"'"e~ — >■ e"'"e~(7) and e"'"e^ r+r^(7) asymmetries. The uncertainty 
on the asymmetry measurement from background subtraction is estimated to be 0.0005. 

The asymmetry is determined from the number of events observed in the forward and 
backward hemispheres, correcting for polar angle dependent efficiencies and background. The 
scattering angle is defined by the polar angle of the electron, 6^-- The determination of the 
asymmetry is repeated defining the angle by the positron, 6^+, and taking the average of the 
two ApB values. This reduces the sensitivity of the result to the size of the interaction region 
and its longitudinal offset. 

Alternatively, we determine the forward-backward asymmetry using the scattering angle in 
the rest system of the final state electron and positron: 

cosr= . (5) 

sm Oe- + sm Oe+ 

This definition minimises the sensitivity to photon emission. A Monte Carlo study shows that it 
differs by less than 0.0005 from the above definition of Ap^ due to different radiative corrections. 
After correcting for this difference in the data the two approaches yield forward-backward 
asymmetries consistent within 0.0015 which is taken as an estimate of the remaining uncertainty 
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of the scattering angle from the knowledge of the interaction point. The contributions of the 
systematic error on the asymmetry measurement are summarised in Table El 

The differential cross sections of the process e"'"e~ — > e"'"e~(7) at three different centre- 
of-mass energy points are shown in Figure |221 together with the prediction of the ALIBABA 
program. 

10 Results on Total Cross Sections and Forward- Backward 
Asymmetries 

The results of the measurements of the total cross section performed between 1993 and 1995 
in the four reactions e"'"e~ — >■ hadrons(7), e"'"e~ — > fi'^ fi~ {'-/) , e+e~ r+r^(7) and e"''e~ — *■ 
e"'"e~(7) are listed in Tables IT^ to 1171 The measured cross sections for e+e~ —>■ hadrons(7) 
are corrected to the full solid angle for acceptance and efficiencies, keeping a lower cut on the 
effective centre-of-mass energy of ^/s' > O.ly^. The measured cross sections for muon and tau 
pairs are extrapolated to the full solid angle and the full phase space using ZFITTER. The 
quoted Bhabha cross sections are for both final state leptons inside the polar angular range 
44° < 6 < 136°, with an acoUinearity angle ^ < 25° and for a minimum energy of 1 GeV of the 
final state fermions. In Table El the s-channel contributions to the cross section extrapolated 
to the full phase space are also given. Their calculation is described in Appendix ^ and they 
can be compared to the measurements of the other leptonic final states fTables IT^ and ITH)). 
Results of the measurements performed between 1990 and 1992 are presented in Reference jl]. 

Figures El to EZI compare the measurements of the total cross sections performed in 1990 — 95 
at the Z pole to the result of the fit to all cross section measurements imposing lepton uni- 
versality described in section 112.11 For Bhabha scattering the contributions from the s- and 
t-channels and their interference are displayed separately. Good agreement between measure- 
ments in different years is observed. 

The measurements of the forward-backward asymmetry performed between 1993 and 1995 
in the leptonic reactions e+e~ — > ;U+/i~(7), e"'"e~ — > r+r^(7) and e+e~ e+e~(7) are listed 
in Tables ITHl to 1201 For muon and tau pairs the results are extrapolated to the full solid angle 
keeping a cut on the acoUinearity of ^ < 15° and ^ < 10°, respectively. The measurements for 
the process e"'"e~ — > e"'"e~ (7) apply to the same polar angular range and cuts as the total cross 
section. Table QUI contains also the s-channel contributions to the asymmetry (see Appendix 
to be compared to the measurements for muon and tau pairs. 

Figures |2H1 to (201 compare these measurements to the results of the fit to all hadronic and 
leptonic cross section and forward-backward asymmetry measurements imposing lepton uni- 
versality. For the Bhabha scattering the difference of the forward and backward cross sections 
in the s- and t-channels and in the interference, all normalised to the total cross section, are 
displayed separately. Good agreement between measurements in different years is observed. 

For the fits presented in the following sections we include the cross section and forward- 
backward asymmetry measurements from 1990—92 |4j. All our measurements at the Z resonance 
performed in the period 1990 — 95 are self-consistent. Qualitatively this can be seen from 
Figure where for all 175 measurements the absolute difference between the measurements 
and the expectations, divided by the statistical error of the measurements, is shown. The 
expected cross sections and forward-backward asymmetries are calculated from the result of the 
five parameter fit presented in Section 112.31 The scattering of our measurements is compared 
with the one expected from a perfect Gaussian distribution. The agreement is satisfactory 
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considering that due to their comphcated correlations, systematic errors cannot be taken into 
account in this comparison. 



11 Fits for Electroweak Parameters 

Different analyses are used to extract electroweak parameters from the measured total cross 
sections and forward-backward asymmetries. 

Firstly, we determine the electroweak parameters making a minimum of assumptions about 
any underlying theory, for example the SM. The first analysis uses only the total cross section 
data to determine the parameters of the Z boson, its mass, the total and partial decay widths 
to fermion pairs. The second analysis also includes the asymmetry data, which allows the 
determination of the coupling constants of the neutral weak current. In a third analysis we fit 
the cross section and forward-backward asymmetry measurements in the S-Matrix ansatz P7j 
where all contributions from 7/Z-interference are determined from the data. Finally, all our 
measurements on electroweak observables are interpreted in the framework of the SM in order 
to determine its free parameters. 



Lowest Order Formulae 

In all analyses, a Breit-Wigner ansatz is used to describe the Z boson. The mass, mz, and 
the total width, Fz, of the Z boson are defined by the functional form of the Breit-Wigner 
denominator, which explicitly takes into account the energy dependence of the total width [38\ . 
The total s-channel cross section to lowest order, cr°, for the process e"'"e~ ff, is given by the 
sum of three terms, the Z exchange, a^, the photon exchange, a°, and the 7/Z-interference, 
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where is the electric charge of the final-state fermion, A*"^ its colour factor, and a the elec- 
tromagnetic coupling constant. The pure photon exchange is determined by QED. 

The first analysis treats the mass and the total and partial widths of the Z boson as free 
and independent parameters. The interference of the Z exchange with the photon exchange 
adds another parameter, the 7/Z-interference term, Jf, besides those corresponding to mass 
and widths of the Z. Since in the SM |cr°^^.(s)| ^ cr°(s) for centre-of-mass energies close to mz, 
it is difficult to measure Jf accurately using data at the Z only. The 7/Z-interference term is 
usually taken from the SM [SHEni, thus making assumptions about the form of the electroweak 
unification. 

The second analysis determines the vector and axial-vector coupling constants of the neutral 
weak current to charged leptons, and (7^, by using the forward-backward asymmetries in 
addition to the total cross sections. In lowest order, for ^/s = mz and neglecting the photon 
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exchange, the s-channel forward-backward asymmetry for the process e"'"e 
by: 



(7) is given 
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The energy dependence of the asymmetry distinguishes gy and EH] • The experimental preci- 
sion on the couphng constants is improved by also including information from tau-polarisation 
measurements which determine Ae and A^- independently. 

In Equation ini the leptonic partial width, F^, and the leptonic 7/Z-interference term, Ji, 
are now expressed in terms of Qy and g^^: 



(4)^ + {Af 



(8) 



where Gp is the Fermi coupling constant. The hadronic cross section is given by the sum over 
the five kinematically allowed flavours and their colour states. Because no separation of quark 
flavours is attempted, this approach cannot be applied to the hadronic final state. Therefore, 
the parameterisation of the first analysis is used to express the hadronic cross section in terms 
of Fhad and Jhad- 

Our data are also interpreted in the framework of the S-Matrix ansatz ^7\, which makes 
a minimum of theoretical assumptions. This ansatz describes the hard scattering process of 
fermion-pair production in e"'"e~-annihilations by the s-channel exchange of two spin-1 bosons, 
a massless photon and a massive Z boson. The lowest-order total cross section, a^^^, and 
forward-backward asymmetry, Ap-^, for e"'"e~ — >■ ff are given as: 
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The S-Matrix parameters rf, jf and gf are real numbers which express the size of the Z 
exchange, 7/Z-interference and photon exchange contributions. Here, rf and jf are treated 
as free parameters while the photon exchange contribution, gf, is fixed to its QED prediction. 
Each final state is thus described by four free parameters: two for cross sections, r\°^ and 
and two for forward-backward asymmetries, r™ and j™. In models with only vector and 
axial-vector couplings of the Z boson, these four S-Matrix parameters are not independent of 
each other: 
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Under the assumption that only vector- and axial- vector couplings exist, the S-Matrix ansatz 
corresponds to the second analysis discussed above without fixing Jhad to the SM. 

The S-Matrix ansatz is defined using a Breit-Wigner denominator with s-independent width 
for the Z resonance. To derive the mass and width of the Z boson for a Breit-Wigner with s- 
dependent width, the following transformations are applied mz = mz + 34.1 MeV and 
Tz = Tz + 0.9 MeV. 



Radiative Corrections 

The QED radiative corrections to the total cross sections and forward-backward asymmetries 
are included by convolution and by the replacement a — > a{s) = a/{l — Aa) to account for 
the running of the electromagnetic coupling constant j^E]. 

Weak radiative corrections are calculated assuming the validity of the SM and as a function 
of the unknown mass of the Higgs boson. The coupling constants which are real to lowest order 
are modified by absorbing weak corrections and become complex quantities |12]. Effective 
couplings, Qy and g^^, are defined which correspond to the real parts. When extracting cjy 
and from the measurements, the small imaginary parts are taken from the SM. Observables 
such as the leptonic partial widths (Equation IH)) and the leptonic pole asymmetry (Equation [Tj) 
are redefined by replacing the vector and axial-vector coupling constants by these effective 
couplings. 

The effective couplings of fermions are expressed in terms of the effective electroweak mixing 
angle, sin^^W) and the effective ratio of the neutral to charged weak current couplings, p = 
1/(1 -Ap) fm-. 

9v = (4 - 2 Qi sin^^w) 

= Vpli, (11) 

where /| is the third component of the weak isospin of the fermion f. Due to weak vertex 
corrections, the definitions of p and sin^6'-w depend on the fermion. However, except for the 
b-quark, these differences are small compared to the experimental precision. Therefore, we 
define sin^^w as the effective weak mixing angle for a massless charged lepton. It is related to 
the on-shell definition of the weak mixing angle, sin^^w by the factor k: 



K = 2n 'With sm fe'w = 1 Y ■ (12) 



Fits in the SM 

The fourth analysis to determine electroweak parameters uses the framework of the SM. By 
comparing its predictions with the set of experimental measurements, it is possible to test the 
consistency of the SM and to constrain the mass of the Higgs boson. 

The input parameters of the SM are a, the fermion masses, mn, mz and the mass of the W 
boson, mw- QCD adds one more parameter, the strong coupling constant, as, which is relevant 
for hadronic final states. The Cabbibo-Kobayashi-Maskawa matrix, relating electroweak and 
mass eigenstates of quarks, is not important for total hadronic cross sections in neutral current 
interactions considered here. Concerning the fermion masses, only the mass of the top quark 
is important for SM calculations performed below. All other masses are too small to play a 
significant role or are known to sufficient precision. 
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Generally, in SM calculations for observables at the Z resonance, the mass of the W is 
replaced by the Fermi coupling constant, Gp, which is measured precisely in muon decay [43\ . 
These two parameters are related by 

Gp ™ 1 1 ^^3^ 



y/2 2 ml sin^e^f^ cos^6w 1 - A 



r 



where Ar takes into account the electroweak radiative corrections. These corrections can be 
split into QED corrections due to the running of the QED coupling constant. Act, and pure 
weak corrections, Ar^ 



Ar = Aa + Ar„ 
Ar„ = - cot^^w Ap + Ar^em ■ (14) 

The corrections Arrem, not absorbed in the p-parameter, are smaller than the main contribu- 
tions discussed below but are nevertheless numerically important [ISllini and included in the 
calculations. 

Weak radiative corrections originate mainly from loop corrections to the W propagator due 
to the large mass splitting in the top-bottom iso-spin doublet and Higgs boson loop corrections 
to the propagators of the heavy gauge bosons jUj. To leading order they depend quadratically 
on rrit and logarithmically on mn. The Zbb vertex receives additional weak radiative corrections 
which depend on the top mass. Through the measurements of weak radiative corrections our 
results at the Z are sensitive to the mass of the top quark and the Higgs boson. This allows to 
test the SM at the one-loop level by comparing the top mass derived from our data with the 
direct measurement and to estimate the mass of the yet undiscovered Higgs boson which is one 
of the fundamental parameters of the SM. With this procedure the relevant parameters in SM 
fits are mz, mt, m^, a{mi) and a^- 



Fitting Programs and Methods 

The programs ZFITTER j48j and TOPAZO ^49| are used to calculate radiative corrections and 
SM predictions. For computational reasons the fits are performed using ZFITTER. 

Both programs include complete 0{a'^) and leading 0{a^\v?{s/m1)) QED calculations of 
initial state radiation [20! • Final state corrections are calculated in 0{a) for QED and 0{a^) [SI] 
for QCD including also mixed terms 0{aas)- Interference of initial and final state radiation is 
included up to 0{a) corrections. Pair production by initial state radiation is implemented j52j. 

Electroweak radiative corrections are complete at the one-loop level and are supplemented 
by leading C(G|mt) and sub-leading 0(G|mtm|) [SB] two-loop corrections. Complete mixed 
QCD-electroweak corrections of 0{aas) with leading 0{G-pTnla1) terms are included to- 
gether with a non-factorizable part [S3] . 

For the reaction e"'"e~ — e"'"e~(7) the contributions from the t-channel photon and Z bo- 
son exchange and the s/t- interference are calculated with the programs ALIBABA [SH] and 
TOPAZO (see Appendix HJ). 

Electroweak parameters are determined in fits using the MINUIT [SZ] program. The 
is constructed from the theoretical expectations, our measurements and their errors, including 
the correlations. Apart from experimental statistical and systematic errors, and theoretical 
errors, we take into account uncertainties on the LEP centre-of-mass energy. Technical details 
of the fit procedure are described in Appendix |Bj 
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Theoretical uncertainties on SM predictions of cross sections, asymmetries, Z decay widths 
and effective couphng constants are studied in detail in Reference [HHlEn]- Errors on the 
theoretical calculations of cross section and forward-backward asymmetries based on total and 
partial Z widths, effective couplings or S-Matrix parameters, as used in the fits of Section ^1 
arise mainly from the finite precision of the QED convolution. They are found to be small 
compared to the experimental precision and do not introduce sizeable uncertainties in the fit 
for Z parameters. Residual SM uncertainties in the imaginary part of the effective couplings 
are even smaller. 

The only exceptions are the theoretical uncertainty on the luminosity determination, dis- 
cussed in Sectional and the treatment of t-channel and s/t-interference contributions to the 
e"'"e~ final state due to missing higher order terms and the precision of the ALIBABA pro- 
gram [nHj. Uncertainties on the Bhabha cross section and forward-backward asymmetry from 
calculations of the t-channel and s/t- interference contributions are discussed in Appendix 1X1 

Additional uncertainties arise in the calculation of SM parameters from the application 
of different re-normalisation schemes, momentum transfer scales for vertex corrections and 
factorisations schemes ^58j. By comparing different calculations as implemented in ZFITTER 
and TOPAZO, we find that the impact of these theoretical uncertainties on the fit results for 
SM parameters presented in Section El is negligible compared to the experimental errors. 

ZFITTER and TOPAZO calculations in the SM framework are performed based on five input 
parameters: the masses of the Z and Higgs bosons, the top quark mass, the strong coupling 
constant at mz and the contribution of the five light quark flavours, Aa||^j, to the running 
of the QED coupling constant to mz- For comparison to the SM we use the following set of 
values and uncertainties [H HIHTHHlj : 

mz = 91189.8 ± 3.1 MeV, mt = 173.8 ± 5.2 GeV , 95.3 GeV < mn < 1 TeV , 
Aaf2d = 0.02804 ±0.00065, = 0.119 ± 0.002 . ^ 

The central values are calculated for mn = 300 GeV. This arbitrary choice is motivated by 
the logarithmic dependence of electroweak observables on mn and it leads to approximately 
symmetric theoretical errors. 

We use the default settings of ZFITTER which provide the most accurate calculations. 
Exceptions are that in all calculations we allow for the variation of the contribution of the five 
light quarks to the running of the QED coupling constant, AaJ^^;^. Secondly, as recommended 
by the authors of ZFITTER, the corrections of Reference j^SI are explicitly calculated for SM 
expectations and in fits in the SM framework (Section IT3|l. In all other cases they are absorbed 
in the definitions of the parameters. 

12 Determination of Z Parameters 
12.1 Mass, Total and Partial Widths of the Z 

We determine the mass, the total width and the partial decay widths of the Z into hadrons, 
electrons, muons and taus in a fit to the measured total cross sections. These parameters 
describe the contribution of the Z exchange to the total cross section. The photon exchange 
and 7/Z-interference contributions are fixed to their SM expectations. Two fits are performed: 
one assuming and one not assuming lepton universality, where in the first one a common 
leptonic width is defined as the decay width of the Z into a pair of massless charged leptons. 
The results of both fits are summarised in Table |^ and the correlation coefficients for the 
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parameters determined in the two fits are given in Tables |22l and ESI respectively. The partial 
decay widths into the three charged lepton species are found to be consistent within errors. It 
should be noted that due to the mass of the tau lepton, is expected to be 0.19 MeV smaller 
than r^. 

Our new results with significantly reduced errors are in agreement with the SM expectations 
and our previous measurements For the mass mz and the total width Fz we obtain: 

mz = 91 189.8 ±3.1 MeV, 

Fz = 2 502.4 ±4.2 MeV. (16) 

These are measurements of the mass of the Z boson with an accuracy of 3.4 • 10~^ and of its 
total decay width of 1.7 ■ 10~^. The contribution to the total errors on mz and Fz from the 
LEP energy is estimated by performing fits to the 1993 — 95 data with and without taking into 
account LEP energy errors. From a quadratic subtraction of the errors of the fitted parameters 
we find Amz(LEP) = 1.8 MeV and AFz(LEP) = 1.3 MeV, in agreement with the estimates 
given in Reference [Oj. 

The impact of the uncertainties on SM parameters on the fit results is negligible. The 
largest effect is an uncertainty on the Z mass of ±0.2 MeV caused by the calculation of the 
7/Z-interference contribution when varying the Higgs and top masses and Aa^^^^j in the ranges 
given in Equation ITKl 

Motivated by the different methods used to obtain the absolute scale of the LEP energy in 
the years 1990 — 92, 1993 — 94 and 1995, resulting in different uncertainties, we determine the 
mass of the Z for these three periods independently. The mass values obtained are consistent 
within their statistical errors. 

To check our results on mz and Fz the fit assuming lepton universality is repeated twice: 
i) using only the leptonic cross sections and ii) using only the e+e~ — >• e"'"e~(7) data. The 
results for the mass and total width obtained this way are mz = 91 198.7 ± 8.2 MeV, Fz = 
2 508.1±13.4 MeV using all three lepton species and mz = 91 177±16 MeV, Fz = 2 497±26 MeV 
when using only Bhabha scattering data. Within the errors, dominated by the statistical errors 
of the measurements, these values are in agreement with those given in Table |^ where the 
e"'"e~ —>■ hadrons(7) cross section measurements contribute most. Also, we conclude that there 
is no significant bias introduced in the determination of the mass and the total width of the Z 
boson by the treatment of the t-channel in Bhabha scattering. 

From the difference of the total width and the partial widths into hadrons and charged 
leptons, including their correlations, the decay width of the Z into invisible particles is derived 
to be 

Finv = 499.1 ±2.9 MeV. (17) 

This number is determined in the fit assuming lepton universality and it is in agreement with our 
direct determination of Finv from cross section measurements of the reaction e"'"e~ — >■ i>i''y{'j) [HS] 
which yields Fi^v = 498 ± 12(stat) ± 12(sys) MeV. 

In the SM, the invisible width is exclusively given by the Z decays into neutrinos and the 
result can be interpreted as the number of neutrino generations N^. Using the SM prediction 
Te/r,y = 0.5021 ± 0.0002 for the ratio of the Z decay width into charged leptons and neutrinos 
we obtain: 

N^ = ^i^) = 2.978 ±0.014. (18) 
This formula is used because the experimental precision on the ratio Finv/F^ is better than that 
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12.2 Limits on Non-Standard Decays of the Z 



From the measurements of total and partial Z decay widths presented in the previous section 
we derive experimental limits on additional Z decay widths not accounted for in the SM. These 
limits take into account experimental and theoretical errors added in quadrature. The latter 
are derived from adding in quadrature the changes in the theoretical predictions when varying 
the SM input parameters by their errors as given in Equation^] This is motivated by the fact 
that these parameters are determined in independent experiments with the exception of the 
mass of the Higgs boson. A value of mn = 1 TeV is used here to calculate the Z widths which 
results in the lowest SM predictions and therefore in conservative limits. 

The 95% confidence level (C.L.) limits on non-standard decay widths, Tg^, are calculated 
using the formula jMj : 
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(19) 



where F'^^p is our experimental result, F^^ the SM expectation for mn = 1 TeV and A the 
combined experimental and theoretical error. 

The limits obtained for the total, hadronic, leptonic and invisible widths, as well as for 
the three lepton species, are summarised in Table |^ Also listed are the differences of our 
measurements and the SM expectations together with their experimental and theoretical 68% 
C.L. errors. It should be noted that the results on the total and partial widths are correlated; 
hence the limits derived in this section cannot be applied simultaneously. 



12.3 Fits to Total Cross Sections and Forward-Backward Asymme- 
tries 

The measured leptonic forward-backward asymmetries are included in the fits. Besides mz and 
Fz the measurements are fitted to the hadronic pole cross section, cTjJg^j, the ratios of hadronic 
to leptonic widths, Ri, and the leptonic pole asymmetries, Ap^, which are defined as: 



U l^TT FcFhad Thad /(O/ ^ , , , , 



The advantage of this parameter set is that the parameters are less correlated than the partial 
widths. Two fits are performed, one with and one without assuming lepton universality. The 
results are listed in Table ESI and the correlation matrices are given in Tables |2ni and UTl 

The 68% C.L. contours in the Ap^ — Ri plane are derived from these fits for the three lepton 
species separately and for all leptons combined (Figure l32|l . In this plot the contour of Ap^ — Rr 
is shifted by the difference in expectation for Rr due to the tau mass to facilitate the comparison 
with the other leptons. Also for the forward-backward asymmetries good agreement among the 
lepton species is observed. Our results are in agreement with the SM expectations. 

From the measurements of the forward-backward pole asymmetries the polarisation param- 
eter, Ai, can be derived for the three individual lepton types as well as the average value. The 
results are listed in Table EHl Because of their relation to the measured pole asymmetry (Equa- 
tion E|) the results for Ae, A^ and Ar are highly correlated. They are compared to Ae and At 
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derived from our measurements of the average and the forward-backward tau-polarisation |67j . 
All measurements are in good agreement and yield an average value of 

= 0.1575 ±0.0067. (21) 

12.4 Vector- and Axial- Vector Coupling Constants of Charged Lep- 
tons 

The effective coupling constants, gy and are obtained from a fit to cross section and forward- 
backward asymmetry measurements, and including our results from tau-polarisation. We use 
the results from tau-polarisation on Ag and Ar as given in Table EHl together with a 8% corre- 
lation of the errors. The inclusion of tau-polarisation results significantly improves the deter- 
mination of the effective coupling constants. 

Fits with and without assuming lepton universality are performed and the vector and axial- 
vector coupling constants so obtained are listed Table 1221 The axial- vector coupling constant of 
the electron is taken to be negative, in agreement with the combination of results from neutrino- 
electron scattering and low energy ApB measurements [HE!- All other signs are unambiguously 
determined by our measurements. 

The 68% CL. contours in the Qy-Qa plane are shown in Figure IHHl revealing good agreement 
among the three lepton species and thus supporting lepton universality in neutral currents. This 
is quantified by calculating the ratio of muon and tau to electron coupling constants, taking 
into account their correlations (see Table EDI)- 

The average vector and axial-vector coupling constants of charged leptons are found to be 

= -0.0397 ±0.0017, = -0.50153 ±0.00053. (22) 

The resulting axial- vector coupling constant g^^ is significantly different from its lowest order SM 
expectation —1/2. This is interpreted as proof for the existence of weak radiative corrections 
from higher order processes and corresponds to a measurement of the p-parameter of 

p = 1.0061 ±0.0021. (23) 

The evidence for the existence of weak radiative corrections is illustrated in Figure EHl 

Measurements of forward-backward asymmetries and tau-polarisation can be compared in 
terms of the effective weak mixing angle, sin^^w, defined by the ratio of the coupling constants 
(Equation HT}. From the average leptonic pole asymmetry and the tau polarisation, the values 
listed in Table |^ are obtained. Our results obtained from the measurement of the forward- 
backward asymmetry of b-quarks [nH] and the measurement of the quark charge asymmetry, 
Qfb, [ZOI are also shown. All four measurements of the weak mixing angle are in agreement 
with each other and the average yields: 

sin^^w = 0.23093 ± 0.00066 . (24) 

12.5 Fits in the S-Matrix Framework 

The programs SMATASY [Tlj together with ZFITTER, ALIBABA and TOPAZO are used for 
the calculation of the theoretical predictions, including QED radiative corrections, of total cross 
sections and forward-backward asymmetries. Further details can be found in Reference jZ21- 
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The results of the fits in the S-Matrix framework with and without imposing lepton uni- 
versahty to the cross sections and forward-backward asymmetries measured at the Z resonance 
are shown in Table ESI The fitted parameters for electrons, muons, taus and hadrons are in 
agreement with each other and with the expectations from the SM. The correlations of the 
parameters as obtained in the two fits are shown in Tables ESI and |33J respectively. 

Large correlations among the parameters are observed. Of particular importance is the 
correlation of —0.95 between the mass of the Z boson and the hadronic interference term, 
It causes an increase of the error on niz with respect to the fits performed in Sections 112.11 
and 112.31 where the 7/Z-interference terms are fixed to their SM expectations. The correlation 
between mz and is illustrated in Figure EH 

Comparing the results on the Z boson mass obtained with the two analyses (Table |^ 
or Eni and Table IH^ good agreement is found. From a quadratic subtraction we estimate 
the additional error on the Z mass arising from the experimental uncertainty on the hadronic 
interference term to be: 

mz = 91 185.2 ± 3.1 ± 9.8(j^°*d) MeV . (25) 

The interference between the photon and the Z is measured to much better precision at 
centre-of-mass energies below or above the Z resonance. By adding our measurements of 
hadronic and leptonic cross sections and forward-backward asymmetries above the Z reso- 
nance this contribution to the error on mz is significantly reduced. This will be reported in a 
forthcoming publication. 



13 Results on SM Parameters 

We interpret our measurements in the framework of the SM to check its consistency by com- 
paring our results to other measurements. The strategy will be to test at first QCD radiative 
corrections in terms of as before we verify weak radiative corrections by comparing the top 
mass derived from our data at the Z resonance to the direct measurement. From our mea- 
surements at the Z, the W mass is determined and compared to our result obtained above the 
W-pair threshold. Finally we use all our measurements of electroweak parameters and include 
the direct measurement of rrit to estimate the mass of the SM Higgs boson. 

Fits are performed to our data to determine the set of SM parameters given in Equation ITHl 
The program ZFITTER is used for SM calculations. In all fits the QED coupling constant 
at the mass of the Z is calculated using the constraint on the contribution from the five light 
quark flavours to the running as obtained in Reference [Mj . 

The input data for the SM fits are our measurements of total cross sections and forward- 
backward asymmetries at the Z resonance performed between 1990 and 1995. In addition, 
our results from tau polarisation (Table EH|) . the effective weak mixing angles from b-quark 
forward-backward and form quark charge asymmetry (Table IHTj). as well as the partial decay 
width into b-quarks -Rb = Tb/rhad = 0.2174 ± 0.0032 [73], are included. 

Firstly, the sensitivity of our data to QCD radiative corrections is exploited to determine 
the strong coupling constant at the mass of the Z boson: 

as = 0.1226 io°:So°60- (26) 

In this fit the result of Section 112.11 is obtained again for mz and Aaj^^^^j remains within the 
imposed constraint. The masses of the top quark and the Higgs boson are free parameters and 
their uncertainties are included in the error on as- 
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The value for as is in good agreement with our determination of the strong couphng constant 
from hadronic event topologies as = 0.1216 ± 0.0017 ± 0.0058 [Zl]. We use this measurement 
as an additional constraint and obtain for the top mass: 

mt = 1971?° GeV. (27) 

This result for the top quark mass is based on our measurements of weak radiative cor- 
rections and their interpretation in the SM framework. The agreement with the direct mt 
measurements by the CDF and DO experiments, rrit = 173.8 ± 5.2 GeV jHI], means that the 
bulk of weak radiative corrections indeed originates from the large mass of the top quark. The 
result for the Higgs boson mass obtained in this fit, logig'^H/GeV = 1.99 Igjg, is in agreement 
with the range allowed by the direct search and the SM (Equation ITHjl. 

From the result of this fit which is based on measurements at the Z resonance a value for 
the mass of the W boson is derived: 

mw = 80.523 ± 0.079 GeV . (28) 

The 68% C.L. contour in the mt-m-w plane obtained in this fit is shown in Figure ESI This 
result for the W mass agrees well with our direct measurements of performed at centre-of- 
mass energies between 161 GeV and 183 GeV [7HH77] which yield a combined value of mw = 
80.61 ±0.15 GeV. 

We include the direct measurement of m-w to determine the electroweak radiative corrections 
Ar, the parameters p and k, the effective weak mixing angle sin^^vv and the on-shell definition 
sin 6'w as well ClS db combined result for mw: 

Ar = 0.0257 ±0.0043, p = 1.0078 ± 0.0017 , k = 1.0493 ± 0.0053 , , . 

sin2^w = 0.23075 ±0.00054, sin^^w = 0.2199 ± 0.0013 , mw = 80.541 ± 0.069 GeV . ^ ' 

This value for the effective weak mixing angle sin^^w derived in a SM fit is in good agreement 
with the result obtained, in a less model-dependent way, from measurements of asymmetries 
(Table EH). 

Finally, we constrain the mass of the top quark to the combined value from the direct 
measurement of DO and CDF. The five SM parameters and their correlations obtained in this 
fit are summarised in Tables EHl and EEl respectively. In particular, for the mass of the yet 
undiscovered SM Higgs boson, we obtain a value and an upper limit: 

mn = 361^9 GeV , 

< 133 GeV 95% C.L. (30) 

Figure EEl shows the 68% and 95% C.L. contours in the mt-mn plane and Figure EZI the 
dependence of the of fit on the Higgs mass from which the upper mass limit is de- 
rived. The result is compatible with the result of our direct search for the SM Higgs boson 
mn > 95.3 GeV 

14 Summary and Conclusion 

We report on the precise measurements of total cross sections and forward-backward asymme- 
tries of the reactions e+e~ — > hadrons(7), e"'"e~ —>■ /i"'"/i~(7), e+e~ —>■ t~^t~{j) and e+e~ — 
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e'^e~(7) at centre-of-mass energies at the peak and the wings of the Z resonance performed in 
the years 1993 — 95. A total luminosity of 103 pb~^ corresponding to 2.5 million hadronic and 
250 thousand leptonic decays of the Z was collected which significantly improve our measure- 
ments of the resonance curve. Including the data samples collected in previous years, the total 
number of Z decays observed by the L3 detector during the first phase of LEP amounts to 4 
million which are used to determine the properties of the Z and other SM parameters. 

All our measurements are consistent with the hypothesis of lepton universality. Prom the 
measured total hadronic and leptonic cross sections we obtain: 

mz = 91 189.8 ± 3.1 MeV , Tz = 2 502.4 ± 4.2 MeV , 
Thad = 1 751.1 ± 3.8 MeV , Te = 84.14 ± 0.17 MeV . 

From these results, the decay width of the Z into invisible particles is derived to be Finv = 
499.1 ± 2.9 MeV, which in the SM corresponds to a number of light neutrino species of: 

A^^ = 2.978 ± 0.014 . (31) 

Including our measurements of leptonic forward-backward asymmetries and tau polarisa- 
tion the effective vector and axial-vector couphng constants of charged leptons to the Z are 
determined to be: 

= -0.0397 ±0.0017, = -0.50153 ±0.00053. (32) 

For the effective weak mixing angle we obtain: 

sin^^w = 0.23093 ± 0.00066 , (33) 

including our measurements of the b-quark forward-backward and quark charge asymmetries. 

Our measurements are sensitive to higher order weak radiative corrections which depend on 
the masses of the top quark and the Higgs boson. Using in addition our measurements of the 
partial width Z — > bb and ttg, we derive in the SM framework a top quark mass 

mt = 1971?^ GeV, (34) 

which is in agreement with the direct measurements of mt. Using our direct measurement of 
mw and the knowledge of mt our data constrain the mass of the Higgs boson to 

ma < 133 GeV 95% C.L. (35) 
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Appendix 



A Treatment of Contributions related to the t-channel 

In the case of the process e^e~ — >■ e^e~ (7) , the existence of the t-channel exchange of photons 
and Z bosons and its interference with the s-channel exchange lead to additional complica- 
tions. Analytical programs to calculate this process, such as ALIBABA and TOPAZO, are not 
directly suited for fitting purposes, as computationally they are very time consuming. Thus, 
the following procedure is adopted. During the initialisation of a fit, ALIBABA is used once 
to calculate the predictions of the t-channel and s/t- interference contributions to the measured 
e"'"e~ — >■ e"'"e~ (7) cross sections and forward-backward asymmetries. Calculations are performed 
at several centre-of-mass energy values in the vicinity of the data points allowing for a reduction 
of the statistical error of the Monte Carlo integration used by ALIBABA and the calculation 
of derivatives needed to construct the covariance matrix (see Appendix El- ZFITTER is em- 
ployed during the fits to calculate the corresponding s-channel contributions as a function of 
the varying electroweak parameters. 

The contribution from s/t-interference depends on the fit parameter, most importantly on 
the Z boson mass. This dependence is taken into account by converting the difference between 
the Z mass used in the initialisation and the current fit value into an equivalent shift in the 
centre-of-mass energy at which the t-channel and s/t-interference contributions are calculated. 
The dependence of the s/t-interference on mz is responsible for a correlation of the results for 
mz and the electron s-channel cross section which amounts to +11% between mz and Re in 
the nine parameter fit (Table EEI)- 

In the analytical program ZFITTER polar angular cuts can only be applied on the positron 
while in the experimental cross section measurements both electron and positron are required 
to lie within the fiducial volume. Correction factors are calculated with TOPAZO which allows 
for both types of cuts. Finally, cross sections of e"'"e~ e+e~(7) to be compared to the 
experimental measurement are calculated as: 

(T*+"/*'^^(44° < ^c-,^c+ < 136°) + a''^^{U° < 6,+ < 136°) R 
a^.T0(44°<g^_^g^^<136°) 

a',TO(AA° < 6^+ < 136°) ^ ^ 

All cross sections are defined for an acoUinearity angle cut of ^ < 25°. The indices AL, TO 
and ZF label cross sections calculated with the ALIBABA, TOPAZO and ZFITTER programs, 
respectively. Cuts on the polar angle of the electron and positron are given in parentheses. This 
procedure is applied because it combines the most accurate treatment of electroweak radiative 
corrections as available in ZFITTER and TOPAZO with the most complete calculations of t- 
channel contributions in ALIBABA. In the case of the forward-backward asymmetry, ApB = 
(cTp — (Tb) / (ctf+ctb), the analogous calculations for the forward, cxp, and backward cross sections, 
(Tb, are performed. 

The cross sections calculated with the ALIBABA and TOPAZO programs for all 1990 — 95 
data sets are listed in in Table EH The theoretical uncertainties on the calculation of the t- 
channel contributions are listed in Table EH These errors are applied to the measured cross 
sections and forward-backward asymmetries in the fits (see Appendix lB|l. 



th 



with R 
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Due to the contribution of the t-channel the results from the Bhabha channel cannot be 
directly compared to the measurements of the other leptonic final states. To permit such a 
comparison, Equation IHEl is used to calculate the s-channel contributions, replacing a^^ by the 
measurements and using ZFITTER for the extrapolation to the full solid angle. The s-channel 
cross sections obtained this way, al, without any cuts, and s-channel asymmetries, A^^, with 
an acoUinearity angle cut of ^ < 25°, are given in Tables IT7I and 1201 



B Construction of the Covariance Matrix 



All fits for electroweak parameters described in this article are performed in the error matrix 
analysis. They consist of minimising a function defined as 



X 



D'^V-^ D. 



(37) 



where D is a column vector with elements defined by the difference between measurements 
Vt^"^ and theoretical expectations Vlf^ which are calculated during the fit as a function of the 
fit parameters: 

Di = nT" - nf" . (38) 

The index i runs over all cross section (f2j = cxj) and forward-backward asymmetry {Qi = Ap^ ^) 
measurements considered in the fit. There are 100 cross section and 75 forward-backward 
asymmetry measurements at the Z resonance taken in 1990 — 95. 

The covariance matrix V is constructed in the following way from all experimental and 
theoretical errors affecting the measurements. The diagonal elements, Vu, are obtained adding 
all individual errors of measurement i in quadrature 



Vu = (/.Ar)' + (/,Ar)V(Ar)'+(Af^)' 



_Tr_lum I T/LEP I T^ecms , yt 



(39) 



where Af^^ is the statistical error of the measurement, A^'^'^ the uncorrelated part of the sys- 
tematic error. A™'' = 6'^°'^ af^ the correlated systematic error which scales with the expected 
value and A^*^*^ the correlated part of the systematic error which does not scale. In case of the 
forward-backward asymmetry the systematic error can simply be expressed in terms of Af^** 
only. The statistical and uncorrelated systematic errors which are derived from the measure- 
ments (Tables ITU to l20|l are scaled during the fit with factors fi to the expected errors using 
the theoretical expectations: 
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(40) 



For cross section measurements there is an additional contribution of the luminosity mea- 
surement, V^^™. Its calculation as well as the contribution from the uncertainty on the LEP 
centre-of-mass energy, V^^^, and its spread, both applied to cross sections and ApB 

measurements, and the theoretical uncertainty on the subtraction of the t-channel and s/t- 
interference contribution to Bhabha scattering, V^*, are described below. 

The off-diagonal elements are constructed from correlated error sources: 
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LEP 



3 

I T / ^cms 



7^ j) 



(41) 
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Experimental systematic errors, A'^™ and A'^'^'^, are only applied to elements connecting the 
same observable, either cross section or asymmetry, of the same reaction. As above, for 
forward-backward asymmetry measurements contributions from A'^°'' and are not ap- 

plicable. Contributions from V^^"™, Vl-^^, V^J"^'' and l^* enter also into off-diagonal elements 
connecting measurements of different observables or reactions. 

All statistical errors and the individual contributions to systematic errors for the measure- 
ments performed in 1993 — 95 are listed in Tables El to EOl For the measurements performed 
in 1990 — 92 systematic errors are quoted in Reference j4j as relative errors, 6i, for the cross 
section and absolute errors. A,, for the forward-backward asymmetries. The correlation among 
these systematic errors is treated by using Vij = {mm{6i,6j)Y af'af^ for the cross sections and 
Vij = (min(Aj, Aj))^ for the asymmetries. Contributions from uncertainties on the luminosity 
and the LEP energy are added, where applicable. 

Correlations between experimental systematic errors in the 1990 — 92 and 1993 — 95 data 
sets are estimated in the same way by using the smaller values of 1993 — 95. Exceptions are 
the measurements of the e+e~ — > t~^t~{'j) cross section where due to the completely revised 
analysis in 1993 — 95 an additional factor of 0.72 is applied to 5™''. Other contributions to 
elements connecting measurements of 1990 — 92 and 1993 — 95 are from V^j"'-^% V^*- and the 
theoretical uncertainty in V^^"™. 

Uncertainty on the Luminosity 

For the 1993 — 95 cross sections the contributions to the covariance matrix from errors on the 
luminosity measurement are obtained from the sum of the total experimental errors, including 
their correlations, and the theoretical uncertainty: 

T/lum / rlum,exp _th\ / rlum,exp th\ lum,exp . / rlum,stat\ ^ Jh th , /rlum,th\^ _th„th /'/lo^ 

Kj = [^k (^i ) [^i (^j ) Pki ' j or^ +[b ' ) a J . (42) 

The indices k and / label the years of the data sets i and j. Total experimental errors on 
the luminosity, (^i^^-'^^p, are given in Tabled] and their correlations, p^kT^''^^^, in Table El The 
statistical error on the luminosity measurement, ^'u^'^^t^t^ is given in Table El for the nine data 
taking periods. It applies only to cross section measurements i and j performed in the same 
period m. 

The combined experimental and theoretical error on the luminosity determination in 1990 — 
92 is 5^"™ = 6%c. It is treated as fully correlated and the corresponding terms in the covariance 
matrix are calculated as Vjl"™ = (5^"™)^ a^^al^. For the 1990 — 92 data the statistical error of the 
luminosity measurement is included in the quoted statistical errors of the hadron cross section 
measurements and it is neglected for the leptonic cross sections. The theoretical uncertainty, 
^ium,th _ Q_gx%(„ is fully correlated among all cross section measurements in 1990 — 95. 

The measurement of the luminosity, and hence of cr^^^, depends on the centre-of-mass energy 
E due to the approximate l/E"^ dependence of the Bhabha cross section: 

J cxp exp 

dE ~ E, ■ ^^^^ 

Because of the 7/Z- interference and higher order contributions the exponent Km differs from 2 
and it is calculated with BHLUMI to be = 1.95, 2.27 and 1.97 for the peak— 2, peak and 
peak+2 data sets, respectively. This dependence causes a small contribution to the uncertainty 
on the cross section measurements and it is taken into account for the 1993 — 95 data together 
with the error on the LEP energy. 
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Uncertainty on the Centre-of-Mass Energy 

The errors on the LEP centre-of-mass energy are transformed into equivalent errors of cross 
section and asymmetry measurements using the partial derivatives of the theoretical cross sec- 
tions and forward-backward asymmetries with respect to the centre-of-mass energy, dQf^/dE. 
The dependence of the measured cross section, via the luminosity, on the centre-of-mass energy 
is included where applicable. 

For the 1993 — 95 data sets the terms in the covariance matrix are determined as 



with 
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(44) 



where k and / are the indices of the LEP energy calibration periods corresponding to the data 
sets i and j, with Hm the appropriate factor defined in Eq. HHl For a pair of data sets taken 
in any of the seven periods with precise LEP energy calibration, vj^f'^ is the corresponding 
element of the LEP energy error matrix as given in Table H For data sets taken during the 
1993 or 1995 pre-scans the LEP energy error is treated as uncorrelated with any other period. 
Hence, only elements connecting data sets taken during the same pre-scan receive contributions 
'^kk^ = (^fc^^)^5 where 5^^^ are the uncertainties on the centre-of-mass energy of 18 MeV and 
10 MeV in 1993 and 1995, respectively. 

The uncertainties on the LEP centre-of-mass energy in 1990 — 92 are listed in TableEl They 
are uncorrelated to the energy errors in 1993 — 95. 

The correction factors for cross sections, f* , and the absolute corrections for forward- 
backward asymmetries, a^cms' cipplied to the measurements to account for the spread of the 
centre-of-mass energy, ecms, can be calculated to sufficient precision from Taylor expansions: 
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(45) 



The contributions to the covariance matrix from the uncertainty on the centre-of-mass energy 
spread, Aecms, are then given by 



T y^cms 

ij 
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A ^cms A £cms 
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'^^T {fl 



2al 



Aer 



for cross sections , 
for ApB ■ 



(46) 



The spread of the centre-of-mass energy and its error for the 1993 — 95 data sets are given 
in Table El The centre-of-mass energy spread used in Reference jl] for the 1990 — 92 data 
has been revised and the values and errors given in Reference jHI are used. The error on the 
centre-of-mass energy spread is fully correlated between all data sets of 1990 — 95, hence all 
elements of the covariance matrix receive a contribution. 
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Error on t-channel Contributions 

The last term in Equations EHl and HU applies to e"*"e~ e"'"e~(7) cross section and Ayb 
measurements only. It accounts for the uncertainty on the calculation of ^-channel and s/t- 
interference contributions: 

V!, = Al A* . (47) 

The indices k and / refer to the centre-of-mass energies of the data set i and j and the corre- 
sponding errors on cross sections and asymmetries, A*, are listed in Table IHHl Off-diagonal 
elements receive a contribution only if the two data sets are both below, on or above the peak. 
Cross section and asymmetry measurements are also connected this way. 

Constraints 

In fits using other results, Q^, than cross section and ApB measurements these additional 
measurements are added with their errors, AJ^, to the function. This applies to our mea- 
surements of tau-polarisation, b-quark and quark charge asymmetries, i?b and as, as well as 
to the value of Aa[^j and the measurement of used in SM fits. Statistical and systematic 
errors are added in quadrature to obtain A^: 

X' = D^V-'D + Y.r\^) . (48) 

The small correlation of 8% between Ae and Ar measured from tau-polarisation is included by 
means of a covariance matrix. 
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1993 




1994 




1995 








peak— 2 


peak 


peak+2 


peak 


peak— 2 


peak 


peak+2 




peak— 2 


12.59 


8.32 


7.45 


5.59 


2.05 


1.80 


1.84 


1993 


peak 




45.69 


7.68 


6.20 


1.69 


1.82 


1.72 




peak+2 






9.57 


5.20 


1.90 


1.96 


2.15 


1994 


peak 








14.30 


1.90 


2.07 


1.92 




peak— 2 










4.49 


2.34 


2.30 


1995 


peak 
peak+2 












30.40 


2.60 
4.15 



Table 1: The covariance matrix, f^; , of the LEP centre-of-mass energy uncertainty 
at the L3 interaction point for the 1994 data set and the two scans of the Z resonance 
performed in 1993 and 1995 jH]. AH values are given in units of MeV^. 





e,ms [MeV] 
1993 1994 1995 


pre-scan 
peak— 2 
peak 
peak+2 


56.8 56.9 
56.6 55.9 

57.0 56.5 56.4 

57.1 56.9 


Aeen,s [MeV] 


±1.1 ±1.1 ±1.3 



Table 2: The spread on the average centre-of-mass energy, Ccms, and its error, Aecms, 
at the L3 interaction point ^ for the nine data sets. The additional scatter from 
the time variation of the mean centre-of-mass energies is included. The uncertainty 
does not depend on the energy. 
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1990 


1991 


1992 


y/s [GeV] 

V L J 


88.231 


89.236 


90.238 


91.230 


92.226 


93.228 


94.223 


91.254 


88.480 


89.470 


90.228 


91.222 


91.967 


92.966 


93.716 


91.294 


88.231 


696 


675 


676 


677 


678 


678 


679 





159 


154 


150 


145 


141 


136 


132 





89.236 




696 


682 


685 


688 


691 


694 





151 


148 


145 


142 


140 


137 


135 





90.238 






706 


692 


698 


703 


709 





142 


141 


141 


140 


139 


139 


138 





91.230 








702 


708 


715 


723 





133 


135 


136 


137 


139 


140 


141 





92.226 










743 


728 


738 





125 


128 


131 


135 


138 


142 


145 





93.228 












764 


753 





116 


122 


126 


133 


137 


143 


148 





94.223 














788 





107 


115 


122 


130 


136 


145 


151 





91.254 
















333 


























88.480 


















93.5 


61.6 


54.1 


44.3 


36.9 


27.1 


19.7 





89.470 




















74.8 


48.2 


40.6 


34.9 


27.3 


21.6 





90.228 






















66.7 


37.8 


33.4 


27.5 


23.0 





91.222 
























45.3 


31.4 


27.7 


24.9 





91.967 


























53.2 


27.9 


26.3 





92.966 




























45.7 


28.3 





93.716 






























57.6 





91.294 






324 



Table 3: The covariance matrix of the LEP centre-of-mass energy uncertainty for 
the 1990 — 92 data sets obtained following references [31j- All values are given in 
units of MeV^. The centre-of-mass energies listed correspond to our measurements 
in Reference j3]. 



Source 


1993 


1994 


1995 


Selection Criteria 


[%c] 


0.48 


0.42 


0.47 


Detector Geometry 


[%c] 


0.63 


0.34 


0.34 


Monte Carlo Statistics 


[%c] 


0.35 


0.35 


0.35 


Total Experimental (6]^^''^^^] 


[%c 




0.86 


0.64 


0.68 


Theory {5,^ ) 


[%c] 


0.61 


0.61 


0.61 


Total Uncertainty 


[%c 




1.05 


0.88 


0.91 



Table 4: Experimental and theoretical contributions to the systematic error on the 
luminosity measurement for different years. Additional contributions to the error 
from statistics and from the uncertainty on the centre-of-mass energy are also taken 
into account in the fitting procedure. 





1993 1994 1995 


1993 
1994 
1995 


1.00 0.59 0.59 
1.00 0.93 
1.00 



Table 5: Correlation coefficients, p^i^''^^^ i between the data sets of different years of 
the total experimental systematic error on the luminosity measurement, as 
given in Table IH 





-^events 




1993 pre-scan 


362 500 


1.66 


peak— 2 


604 535 


1.29 


peak 


651931 


1.24 


peak+2 


588 962 


1.30 


1994 peak 


3129 424 


0.57 


1995 pre-scan 


480 342 


1.44 


peak— 2 


541580 


1.36 


peak 


283 887 


1.88 


peak+2 


554 371 


1.34 


peak combined 


4 908 084 


0.45 



Table 6: Number of events used for the measurement of the total luminosity in the 
nine data taking periods, iVevents, and the corresponding contributions to the error 
of cross section measurements, (5^"™''^*^'^. The numbers correspond to the luminosity 
used for the hadron cross section measurements and are used for leptonic channels 
too. The last line is the sum of the five data sets taken at the peak, indicating this 
error contribution to the measurements of the pole cross sections. 
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Source 


lyyo — 


lyyo 


Monte Carlo statistics 


[%c\ 




n 1 n 




W 1 

[%c\ 




n 91 

U.Zl 


QplppfiriTi pnf«3 


[/ccj 




u.ou 


Trigger 


[%o] 




0.12 


Total scale 


[%c] 


0.39 - 


0.40 


Non-resonant background [pb] 






Detector noise 


[pb] 




1 


Total absolute 


[pb] 


3.2 



Table 7: Contributions to the systematic uncertainty on the cross section e+e~ 
hadrons(7). Except for the contribution from Monte Carlo statistics all errors are 
fully correlated among the data sets. The resulting correlated scale error is 5'^°'' = 
0.39%c. 



Source 


1993 


1994 


1995 


Monte Carlo statistics [%c] 


0.9 


- 1.5 


0.4 


1.7 


- 2.4 


Acceptance 


[%o] 




2.7 


2.7 




3.2 


Selection cuts 


[%c] 




1.3 


1.3 


1.4 


-2.2 


Trigger 


[%c\ 




0.6 


0.6 


0.5 


-0.7 


Resonant background [%c] 




0.3 


0.3 




0.3 


Total scale 


[%o] 


3.2 


-3.4 


3.1 


3.9 


-4.6 


e+e~ e+e" /x"*"//" 


[pb] 










0.1 


Cosmic rays 


[pb] 




0.3 


0.3 




0.3 


Total absolute 


[pb] 


0.3 


0.3 


0.3 



Table 8: Contributions to the systematic uncertainty on the cross section e+e~ 
H'^ IJ,~ {'f) . Except for the contribution from Monte Carlo statistics, all errors are 
fully correlated among the data sets yielding a correlated scale error of 5'^™' = 3.1%c 
for 1993 — 94 data. For the 1995 data this error is estimated to be 3.6%c and it is 
taken to be fully correlated with the other years. 



Source 


1993 - 94 


1995 


Fit procedure 


< 0.0003 


< 0.0003 


Detector asymmetry 


0.0006 


0.0006 


Charge confusion 


0.0004 


0.0010 


Momentum reconstruction 


0.0004 


0.0009 


Background 


0.0001 - 0.0005 


< 0.0003 


Total uncertainty 


0.0008 - 0.0009 


0.0015 



Table 9: Contributions to the systematic uncertainty on the forward-backward asym- 
metry in e"'"e~ — > /i"'"/x~(7). The total uncertainty is assumed to be fully correlated 
among the data sets. 
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Source 


1993 


1994 


1995 


Monte Carlo statistics 


[%c] 


2.3 


-4.2 


1.1 


1.5 


- 1.7 


Tau branching fractions 


[%c] 




2.0 


2.0 




2.0 


Selection cuts 


[%c] 


5.3 


-8.0 


6.0 


6.4 


- 7.5 


Trigger 


[%o] 




0.2 


0.2 




0.2 


Resonant background 


[%c] 


1.4 


-3.3 


1.0 


1.3 


-3.0 


Total scale 


k/ 1 

/f M 


(3.8 


- 9.5 


6.5 


7.5 


- 8.0 


e+e~ — > e+e" r+r~, e+e" 


hadrons [pb] 






0.6 






Cosmic rays 


[pb] 






1.0 






Total absolute 


[pb] 


1.2 



Table 10: Contributions to the systematic uncertainty on the cross section e+e^ 
T+T~(7). The total absolute and 5.7%c of the total scale error are assumed to be 
fully correlated among all data sets. 



Source 


1993 - 1995 


Detector asymmetry 
Charge confusion 
Background 
Helicity bias 


0.0030 
< 0.0001 
0.0010 
0.0004 


Total uncertainty 


0.0032 



Table 11: Contributions to the systematic uncertainty on the forward-backward 
asymmetry in e+e^ r+r~(7). The total uncertainty is assumed to be fully corre- 
lated among the data sets. 



Source 


1993 


1994 


1995 


Monte Carlo statistics 


%o] 


0.4 


- 1.0 


0.4 


0.4 


* Generator 


%o] 




1.0 


1.0 


1.0 


* e-7 discrimination 


%c\ 




0.6 


0.6 


1.5 


* z-vertex position 


%c] 




0.5 


0.5 




* BGO geometry 


%c] 




0.5 


0.5 




* fiducial volume 


%o] 








1.5 


Background 


%o] 


0.7 


- 1.0 


0.4 


0.4 


Selection cuts 


%o] 


1.6 


- 2.7 


0.8 


1.5 


Total scale 


%0 




2.3 


-3.2 


1.7 


2.8 



Table 12: Contributions to the systematic uncertainty on the cross section e"'"e — * 
e^e~(7). Error sources indicated as * are correlated yielding a total correlated scale 
error of 5^°' = lA%c for the 1993 - 94 data and 2.3%c for the 1995 data. 
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Source 


1993 - 1995 


Charge confusion 

2;-vertex 

Background 


0.0020 
0.0015 
0.0005 


Total uncertainty 


0.0025 



Table 13: Contributions to the systematic uncertainty on the forward-backward 
asymmetry in e"'"e~ — >■ e"'"e~(7). 
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x/s [GeVl 


■^^G VG ntr S 


^ [ir J 


a [nb] 


A;^"'^ [nbl 


91.3217 


158 736 


5.21 


30.665 ±0.077 


0.003 


89.4498 


83 681 


8.32 


10.087 ±0.035 


0.001 


91.2057 


281359 


9.34 


30.309 ±0.057 


0.003 


93.0352 


121926 


8.79 


13.909 ±0.040 


0.001 


1993 Totals 


645 702 


31.66 




91.2202 


1359 490 


44.84 


30.513 ±0.026 


0.001 


91.3093 


209 195 


6.90 


30.512 ±0.066 


0.003 


89.4517 


75102 


7.46 


10.081 ±0.037 


0.001 


91.2958 


123 791 


4.08 


30.493 ±0.086 


0.003 


92.9827 


117555 


8.28 


14.232 ±0.041 


0.001 


1995 Totals 


525 643 


26.72 




Total sum 


2 530 835 


103.21 





Table 14: Average centre-of-mass energies, number of selected events, integrated lu- 
minosities and measured cross sections with statistical errors for e"'"e~ — >■ liadrons(7). 
The cross sections are quoted for x/s' > O.ly/s. Apart from the uncorrelated part 
listed, A"'^'^, systematic errors consist in addition of a fully correlated multiplicative 
contribution, 6^°^ = 0.39%c and an absolute uncertainty, A^*^*^ = 3.2 pb. Systematic 
errors from the luminosity measurement (Tables |3] and ^ have to be added. 



[GeV] 


A^events 


£ [pb-i] 


a [nb] 


A,"°^ [nb] 


91.3216 


5134 


5.00 


1.504±0.021 


0.002 


89.4499 


2 739 


8.03 


0.504 ±0.010 


0.001 


91.2056 


9 237 


9.15 


1.482 ±0.016 


0.001 


93.0357 


3 679 


7.97 


0.684±0.011 


0.001 


1993 Totals 


20 789 


30.15 




91.2197 


41 768 


41.38 


1.484 ±0.007 


0.001 


91.3090 


5 772 


5.09 


1.467±0.020 


0.003 


89.4517 


2 789 


7.36 


0.490 ±0.009 


0.001 


91.2949 


3 967 


3.47 


1.483 ±0.024 


0.003 


92.9825 


4 305 


7.83 


0.703 ±0.011 


0.002 


1995 Totals 


16 833 


23.75 




Total sum 


79 390 


95.29 





Table 15: Same as Table IT^ for e+e^ /i+/i~(7). The cross sections are extrapo- 
lated to the full phase space. Apart from the uncorrelated part listed, systematic 
errors consist in addition of a fully correlated multiplicative contribution, 6^°^ = 3.1%o 
(3.6%c for 1995 data), and an absolute uncertainty, A^*^*^ = 0.3 pb. 
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[GeV] 


-'^events 


L [pb-i] 


o [nb] 


A^" [nb] 


91.3221 


4 805 


5.08 


1.492 ±0.021 


0.006 


89.4500 


2 706 


8.08 


0.509 ±0.010 


0.004 


91.2058 


8 506 


9.00 


1.474 ±0.017 


0.007 


93.0358 


3 637 


7.93 


0.718 ±0.012 


0.003 


1993 Totals 


19 654 


30.09 




91.2197 


41439 


43.53 


1.472 ±0.007 


0.005 


91.3096 


7314 


7.75 


1.474 ±0.017 


0.007 


89.4518 


2 352 


7.48 


0.483 ±0.010 


0.002 


91.2951 


3 509 


3.67 


1.503 ±0.025 


0.009 


92.9828 


3 723 


8.25 


0.707±0.012 


0.004 


1995 Totals 


16 898 


27.15 




Total sum 


77991 


100.77 





Table 16: Same as Table ITHfor e"'"e~ r+r^(7). The cross sections are extrapolated 
to the full phase space. Apart from the uncorrelated part listed, systematic errors 
consist in addition of a fully correlated multiplicative contribution, bf^ = 5.7%c and 
an absolute uncertainty, A^'^*' = 1.2 pb. 



[GeV] 


-^events 


C [pb-i] 


a [nb] 


A^'^^ [nb] 


(T^ [nb] 


91.3213 


5 267 


5.19 


1.017±0.014 


0.003 


1.438 ±0.023 


89.4497 


4610 


8.28 


0.566 ±0.008 


0.001 


0.515 ±0.014 


91.2057 


9 834 


9.23 


1.075 ±0.011 


0.002 


1.497 ±0.018 


93.0358 


3 610 


8.39 


0.431 ±0.007 


0.001 


0.703 ±0.012 


1993 Totals 


23 321 


31.09 






91.2197 


43 300 


40.64 


1.075 ±0.005 


0.001 


1.501 ±0.008 


91.3106 


8 200 


7.88 


1.042 ±0.012 


0.002 


1.476 ±0.019 


89.4517 


3 891 


6.91 


0.564 ±0.009 


0.001 


0.511±0.015 


91.2960 


4 310 


4.02 


1.072 ±0.017 


0.002 


1.520 ±0.027 


92.9828 


3 405 


7.90 


0.427 ±0.007 


0.001 


0.696 ±0.012 


1995 Totals 


19 806 


26.72 






Total sum 


86 427 


98.45 







Table 17: Average centre-of-mass energies, number of selected events, integrated 
luminosity and measured cross sections with their statistical and uncorrelated sys- 
tematic errors in the fiducial volume 44° < 9 < 136° and for C, < 25° for the reaction 
e"'"e~ — > e"'"e~(7). The systematic error consists in addition of a fully correlated 
multiplicative contribution, 6f°^ = lA%c (2.3%c in 1995). In the rightmost column 
the s-channel contribution to the total cross section in the full solid angle is listed 
with the statistical error (see Appendix Ej). 
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N 

^ ' events 


^FB 


91.3217 


5 385 


0.009 ±0.015 


89.4497 


2 631 


-0.182 ±0.020 


91.2054 


9150 


0.000 ±0.011 


93.0352 


3 635 


0.119±0.017 


91.2203 


43 416 


0086 ±0 0051 


91.3090 


5 643 


0.022 ±0.013 


89.4517 


2 686 


-0.175 ±0.019 


91.2949 


3 858 


0.030 ±0.016 


92.9825 


4193 


0.104 ±0.015 



Table 18: Results on the forward-backward asymmetry, including an acoUinearity 
cut ^ < 15°, for e"'"e~ — *• /i+/i~(7) together with the average centre-of-mass energies 
and the number of events. The errors are statistical only and a correlated absolute 
systematic error of Af^' = 0.0008 (0.0015 for 1995 data) has to be added. 



[GeV] 


-^events 


^FB 


91.3221 


3 655 


-0.003 ±0.017 


89.4500 


2 090 


-0.138 ±0.022 


91.2058 


6 669 


0.020 ±0.013 


93.0358 


2 822 


0.133 ±0.019 


91.2202 


36 509 


0.0062 ±0.0053 


91.3096 


6317 


0.044 ±0.013 


89.4518 


2 020 


-0.134 ±0.023 


91.2951 


3017 


0.001 ±0.018 


92.9828 


3 263 


0.134 ±0.017 



Table 19: Same as Table fTHl for e"'"e t^t (7). The cut on the acoUinearity angle 
is ^ < 10°. The errors are statistical only and a correlated absolute systematic error 
of Af " = 0.0032 has to be added. 
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x/i [GeV] 


-^events 


^FB 


^FB 


91.3213 


4 009 


0.083 ±0.016 


-0.027 ±0.022 


89.4497 


3 434 


0.311 ±0.016 


-0.119 ±0.039 


91.2057 


7330 


0.111±0.012 


-0.003 ±0.016 


93.0358 


2 679 


0.101 ±0.019 


0.110 ±0.024 


91.2197 


31636 


0.1213 ±0.0056 


0.014 ±0.008 


91.3106 


7861 


0.101±0.011 


-0.001 ±0.016 


89.4517 


3 722 


0.257 ±0.016 


-0.244 ±0.039 


91.2960 


4 083 


0.080 ±0.016 


-0.030 ±0.023 


92.9828 


3 203 


0.055 ±0.018 


0.054 ±0.022 



Table 20: Measured forward-backward asymmetry in the fiducial volume 44° < 9 < 
136° and for ^ < 25° for the reaction e^e~ e^e~(7). The errors are statistical only 
and a correlated absolute systematic error of Af = 0.0025 has to be added. In the 
rightmost column the s-channel contribution to the forward-backward asymmetry 
in the full solid angle is listed for C, < 25° with the statistical error (see Appendix E}- 
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Parameter 


Treatment of Charged Leptons 


Standard 




non-universality 


universality 


Model 


mz [MeV] 


91189. 7±3.1 


91 189.8±3.1 


— 


Tz [MeV] 


2 502.4 ±4.2 


2 502.4 ±4.2 


2 492.7 til 


rhad[MeV] 


1750.9 ±4.7 


1751.1±3.8 


1 739.8 tli 


Te [MeV] 


84.16 ±0.22 


— 


83.91 tg-L^ 


T^. [MeV] 


83.95 ±0.44 




83.91 to°-i4° 


[MeV] 


84.23 ±0.58 




83.72 1°-}4° 


r, [MeV] 




84.14±0.17 


83.91 t|ii4° 


xVdof 


91/94 


91/96 





Table 21: Results of the fits to the 1990 — 95 total cross section data with and without assuming 
lepton universality. 





mz 


Tz 


Thad 




r. 




mz 


1.00 


0.06 


0.14 


-0.03 


0.07 


0.05 


Tz 




1.00 


0.54 


0.55 


0.28 


0.21 


Thad 






1.00 


-0.29 


0.48 


0.36 










1.00 


-0.14 


-0.11 












1.00 


0.19 














1.00 



Table 22: Correlation matrix for the six parameter fit without assuming lepton 
universality in Table 1211 





mz 


Tz 


Thad 




mz 


1.00 


0.07 


0.12 


0.03 


Tz 




1.00 


0.68 


0.71 


Thad 






1.00 


0.12 










1.00 



Table 23: Correlation matrix for the four parameter fit assuming lepton universality 
in Table EH 
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pexp _ pSM pgY] 


[MeV] 




14.6 ±4.2 ±1.7 


22.0 


Thad 


15.2 ±3.8 ±1.5 


21.9 


r 


-1.7±2.9 ±0.23 


4.8 




0.37±0.17±0.05 


0.66 




0.38 ±0.22 ±0.05 


0.75 




0.17 ±0.44 ±0.05 


0.99 




0.64 ±0.58 ±0.05 


1.64 



Table 24: The one-sided upper limits (95% CL.) on non-SM contributions to the 
Z widths, Tg^, as derived from SM calculations and our measurements. Also given 
are the differences of our measurements and the SM expectations for mn = 1 TeV, 
rexp_pSM^ with the experimental and theoretical errors. The results on the total and 
partial widths are correlated and hence the limits cannot be applied simultaneously. 
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Parameter 


Treatment of Charged Leptons 


Standard 




non- universality 


universality 


Model 


mz [MeV] 


91 189.8±3.1 


91 189.5±3.1 




Tz [MeV] 


2 502.5 ±4.2 


2 502.5 ±4.2 


2 492.7^5:2 


(^Ld [nb] 


41.535 ±0.055 


41.535 ±0.055 


41.476 ±0.012 


Re 


20.816 ±0.089 




20.733 ± 0.018 




20.861 ±0.097 




20.733 ± 0.018 


Rt 


20.792 ±0.133 




20.780 ± 0.018 


Re 




20.810 ±0.060 


20.733 ± 0.018 


^FB 


0.0106 ± 0.0058 




0.0151 ± 0.0012 


^FB 


0.0188 ±0.0033 




r\r\ir'i i r\ r\r\-\ c\ 

0.0151 ± 0.0012 


aO,t 
^FB 


0.0260 ±0.0047 




0.0151 ±0.0012 


^FB 




0.0192 ±0.0024 


0.0151 ±0.0012 




158/166 


163/170 





Table 25: Results on the mass, total width, the hadronic pole cross section, crJI^^, 
the ratios of hadronic to leptonic widths, i?^, and the leptonic pole asymmetries, 
ApB, determined from cross section and forward-backward asymmetry data with and 
without assuming lepton universality. The SM expectations are calculated using the 
parameters listed in Equation 1151 





mz 


Tz 


'^had 


Rc 




Rt 


^FB 


^FB 


4U,r 
^FB 




1.00 


0.07 


-0.01 


0.11 


0.00 


0.00 


-0.05 


0.05 


0.03 


Tz 




1.00 


-0.35 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 








1.00 


0.07 


0.08 


0.05 


0.01 


0.01 


0.00 


Re 








1.00 


0.03 


0.02 


-0.15 


0.02 


0.01 


Rfi 










1.00 


0.02 


0.00 


0.00 


0.00 


Rt 












1.00 


0.00 


0.00 


0.01 


^0,6 

^FB 














1.00 


0.01 


-0.01 


AO,^l 

^FB 
















1.00 


0.01 


^FB 


















1.00 



Table 26: Correlation matrix for the nine parameter fit not assuming lepton univer- 
sality in Table 123 
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nil 


J- z 


"^had 




.6/ 
^FB 


777.Z 


1 nn 


n Dfi 


n ni 

U.Ul 


n n? 

U.U 1 


n OA 

u.u^ 


J- Z 




1 nn 

l.UU 


— n 


n nn 

u.uu 


n n9 


'^had 






1 nn 

i .uu 


n 1 9 

U. IZ 


n ni 

U.U i 


Ri 








1.00 


-0.02 












1.00 



Table 27: Correlation matrix for the five parameter fit assuming lepton universality 
in Table 123 



Forward-backward asymmetry 


A, 

Ar 


r, 1 1 Q +0.029 

U-iiy -0.039 

0.210 

0.2911°:^^^ 


Ae 


0.160 ±0.010 


Tau-polarisation 


Ae 


0.1678 ±0.0130 


Ar 


0.1476 ±0.0108 


Average 




0.163 ±0.011 


A, 


0.153 ±0.029 


Ar 


0.152 ±0.010 


Ae 


0.1575 ±0.0067 



Table 28: Measurements of the polarisation parameter A£ from forward-backward 
asymmetries for the three leptons and the combined value. Also listed are the 
results for Ae and Ar obtained from our measurements of tau-polarisation (nZj. At 
the bottom the averages from the combined fit are given. 



Parameter 


Treatment of Charged Leptons 


Standard 




non- universality 


universality 


Model 




-0.0412 ±0.0027 






-g'v 


-0.0386 ±0.0073 






9v 


-0.0384 ±0.0026 






fv 




-0.0397± 0.0017 


0.0358 ±0.0014 


9l 


-0.5015 ±0.0007 






9l 


-0.5009 ±0.0014 






9l 


-0.5023 ±0.0017 






9i 




-0.50153 ±0.00053 


-0 5ni 1 3 +0 00034 
U.OUiit) -0.00022 



Table 29: Results for the vector and axial- vector coupling constants of charged lep- 
tons as obtained from fits to the total cross sections and forward-backward asym- 
metries and including our results from tau polarisation. The errors are correlated. 
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9y/9y 
9y/9y 


0.94 ±0.21 
0.93 ±0.09 


9 a/ 9a 
9 a/ 9 a 


0.9988 ±0.0033 
1.0017±0.0038 



Table 30: Ratios of vector and axial-vector coupling constants obtained from a fit 
to the total cross sections and forward-backward asymmetries including our results 
from tau polarisation. 



Input data 


sin^6'w 


^FB 

Tau polarisation 

.0,b 
^FB 

Qfb 


0.2299 ±0.0013 
0.2304 ±0.0011 
0.2318 ±0.0013 
0.2327 ±0.0017 


Average 


0.23093 ±0.00066 



Table 31: Determination of the effective weak mixing angle sin^^w from different 
reactions: the leptonic forward-backward pole asymmetry, Ap^, tau polarisation 
f Table 1^ . b-quark pole asymmetry [69J, A^^, and the quark charge asymmetry [70] . 
Qfb- The SM prediction is sin^w = 0.23215 
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Parameter 


Treatment of Charged Leptons 


Standard 




non-universality 


universality 


Model 


mz [MeV] 


91 189.3 ±11.2 


91 185.2 ±10.3 


— 


Tz [MeV] 


2502.8±4.6 


2 503.1 ±4.5 


2 492.7111 


^tot 
'had 


2.986 ±0.010 


2.986 ±0.010 


Ql^8 +U.UU9 
Z.yOO _o.012 


».tot 
e 


0.14316 ±0.00089 


— 




».tot 


0.14302 ±0.00082 


— 




^tot 

T 


0.14386 ±0.00104 


— 






— 


0.14336 ±0.00066 


U.i^Z^O -0.00049 


■tot 
Jhad 


0.21±0.63 


0.44 ±0.59 


n 91 +0.0086 
U.ZlOO _o.0093 


■tot 
Jc 


-0.029 ±0.054 







■tot 

JfJ. 


0.035 ±0.046 






^tot 
Jt 


0.073 ±0.048 






^ Lot 




0.045 ±0.035 


0.00409 ± 0.00032 


' e 


0.00174 ±0.00113 








0.00341 ±0.00066 






„FB 
r 


0.00456 ±0.00093 






„FB 




0.00333 ±0.00048 


0.00255 ± 0.00023 


Je 


0.698 ±0.080 






aFB 
Jfi 


0.820 ±0.047 






jFB 
Jt 


0.754 ±0.055 






aFB 

Je 




0.777 ±0.033 


0.7986 ;g:oT2^ 




153/159 


161/167 





Table 32: Results of the fits to the data taken at the Z resonance within the S-Matrix 
framework with and without the assumption of lepton universality. 
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Tz 


^tot 
had 


tot 
c 


^tot 


^tot 
r 


^•tot 
•^had 


•tot 
Je 


^•tot 

Jfl 


•tot 

J T 


FB 

e 


FB 


FB 

T 


Je 


Jfl 


Jt 




1.00 


-0.38 


-0.39 


-0.53 


-0.23 


-0.18 


-0.96 


-0.53 


-0.63 


-0.60 


-0.20 


0.22 


0.14 


0.01 


-0.01 


0.01 


Zj 




1.00 


0.93 


0.64 


0.56 


0.44 


0.41 


0.20 


0.27 


0.26 


0.08 


-0.07 


-0.04 


-0.01 


0.04 


0.03 


tot 
' had 






1.00 


0.66 


0.57 


0.45 


0.43 


0.20 


0.28 


0.27 


0.09 


-0.07 


-0.04 


-0.01 


0.04 


0.03 


tot 
' e 








1.00 


0.40 


0.31 


0.54 


0.32 


0.36 


0.34 


0.21 


-0.12 


-0.07 


-0.03 


0.03 


0.01 


^tot 
M 










1.00 


0.27 


0.24 


0.11 


0.25 


0.16 


0.05 


-0.02 


-0.02 


-0.01 


0.08 


0.02 


^tot 












1.00 


0.19 


0.09 


0.13 


0.20 


0.04 


-0.03 


0.01 


-0.01 


0.02 


0.09 


•tot 
Jhad 














1.00 


0.53 


0.63 


0.60 


0.20 


-0.21 


-0.14 


-0.01 


0.01 


-0.01 


•tot 
Je 
















1.00 


0.35 


0.33 


0.10 


-0.12 


-0.08 


0.17 


0.00 


0.00 




















1.00 


0.39 


0.13 


-0.11 


-0.09 


-0.01 


-0.04 


0.00 


•tot 
Jt 




















1.00 


0.12 


-0.13 


-0.07 


-0.01 


0.00 


-0.04 


^FB 
' e 






















1.00 


-0.05 


-0.03 


0.03 


0.00 


0.00 


».FB 
























1.00 


0.04 


0.00 


0.17 


0.00 


^FB 


























1.00 


0.00 


0.00 


0.15 


•FB 

J ^ 




























1.00 


0.00 
1.00 


0.00 
0.00 
1.00 



Table 33: Correlation of the S-Matrix parameters listed in Table IH^ not assuming 
lepton universality. 







J- / 


tot 
' had 


tot 


J had 


■tot 

Jl 




Ji 


11 ill 


1 nn 




— n 35 


— n 41 


— n Q5 


— n 76 

U . 1 u 


n 14 


— n n3 

u . uo 


-L Z 




1 nn 


n Q3 


n 74 


n S8 


n 9Q 


— n ns 

u . uo 


n nfi 

U . UU 


^tot 
'had 






1 nn 


n 7fi 


n 3Q 


n sn 

u .ou 


— n ns 

u. uo 


n nfi 

u.uu 










1 on 


U.HiO 


u .oo 


— n m 

U. U X 


n OR 

u .uo 


■tot 
Jhad 










i.nn 


n.76 


-n.i4 


n.ns 


■tot 
Jl 












I.nn 


-n.ng 


n.n2 


_,FB 














I.nn 


n.i4 


3i 
















I.nn 



Table 34: Correlation of the S-Matrix parameters listed in Table IH^ assuming lepton 
universality. 



SM parameter 


Fit result 


mz [MeV] 
mt [GeV] 
logio "^h/ GeV 


91 189.n±3.1 
177.4 ±4.8 
1.56±n.33 

n.i226±n.nn4n 
n.n2787±n.nnn62 



Table 35: Results of the fit for SM parameters. In addition to our measurements, 
constraints on mt and Aa[j^^^ as given in Equation HSl are used. 





mz 


mt 


logio "^n/GeV 








I.nn 


-n.ni 


-n.n2 


n.ni 


n.n2 


mt 




I.nn 


n.n8 


-n.n7 


n.ii 


logio mn/G(N 






I.nn 


-n.n3 


-n.i4 










I.nn 


n.n8 












I.nn 



Table 36: Correlation matrix for the fit for SM parameters in Table l35l 
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a/s [GeVJ 


^t+s/t, Ah r^Kl 
CT ' ' [pbj 


K 


t+s/t,AL r 1 1 




-Kb 


lyyu 




OO^^ r 

zzo.o 


n nno 
0.99o5 


1 00 

loo. 2 


n nnoo 
U.9922 


n nn /1 
0.994O 




on oo/^ 


n V 
2oU. 1 


n nn yi 
U.994D 


ono 7 
2U2. 1 


n nno a 
U.99o4 


n nn 
0.995O 




9U.ZO / 


O^^ 1 

zol.o 


n nn 
0.9955 


oil 7 
211. 1 


n nn /1 
U.9940 


n nn^^ 
0.99o5 




9I.Z0U 


1 /I 


n nn ^ 
0.995 / 


1 1 
II0.5 


n nn /i n 
U.9949 


n nn^^ 
0.99o5 




no OOP 


1 ^ ^ 
17.7 


n nn /i 

0.994-3 


1 /I c 
14.5 


n nn /I a 

U.9944 


n nn /I 

0.9940 




yo.zzo 


0.0 


n QQ1 1 
u.yyii 


/I 

0.4 


n QQ99 

u.yyzz 


u.yoyo 




n/1 99Q 

y4. zzo 


1 v 




1 /I c; 
14.0 


n nfifi7 


n nOOfi 

U.9ooo 


1991 


m c 

91. ZOO 


1 /1 
14Z.Z 


n nn c ^ 

11.9957 


1 1 /I 
114.8 
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Table 37: Contribution from the t-channel and s/t-interference and corrections to the 
process e"'"e~ e"'"e~ (7) as calculated with the ALIBABA and TOPAZO programs 
for all centre-of-mass energies: a^+^f^'-^'^ is the contribution to the total cross section 
and R the correction factor for the different polar angle cuts (see Appendix^. The 
corresponding differences of forward and backward cross section, a^^^^'^^ , and the 
correction factors for the forward and backward parts of the cross section, Rp and 
Rb, are also given. 



[GeV] 


88.45 89.45 90.20 91.19 91.30 91.95 93.00 93.70 


Cross section [pb] 
Asymmetry 


1.1 1.2 1.2 1.0 1.0 1.2 1.0 1.0 
0.0012 0.0014 0.0007 0.0007 0.0006 0.0001 0.0002 0.0002 



Table 38: Uncertainties of the ALIBABA program on the calculation of the t- 
channels and s/t-interference contribution to e"'"e~ — >■ e"'"e~(7) To account 

for the smaller fiducial volume used in our analysis the errors are scaled by a factor 
of 0.8. 
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Figure 1: The distributions of the energies measured in the luminosity detectors for small angle 
Bhabha candidates in 1993. The top plot contains the most energetic cluster, £^max, and the 
lower plot shows the energy of the cluster on the opposite side. All selection cuts except the 
one under study are applied. In this and the following figures, the dots are the data and the 
histograms represent Monte Carlo simulations. The vertical arrows indicate the positions of 
the selection cuts (see text). r.-. 
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Figure 2: The distribution of the coplanarity angle A0 for small angle Bhabha candidates. 
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Figure 3: The polar angle distribution of small angle Bhabha events used for luminosity mea- 
surement as observed in the two detectors at —z and +2;. The structure seen in the central 
part of the +2; side is due to the flare in the beam pipe on this side. 
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Figure 4: The distribution of the photon energy, E^, as measured in the luminosity monitors, 
normahzed to the beam energy, £"5. 
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Figure 5: The distribution of the visible energy normahzed to the centre-of-mass energy for 
e+e~ — > hadrons(7) candidates collected in 1994. In this and the following figures the data 
are presented as dots, the Monte Carlo simulations of the signal as open and of the different 
background sources as shaded histograms. 




Figure 6: The distribution of the longitudinal energy imbalance for e"'"e — > hadrons(7) candi- 
dates. 
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Figure 7: The distribution of the transverse energy imbalance for e"'"e — > hadrons(7) candi- 
dates. 
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Figure 8: The distribution of the number of energy clusters in the calorimeters for e 
hadrons(7) candidates in the barrel region (| cos^t| < 0.74). 
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Figure 9: Same as Figure |H1 for events in the endcap region (| cos6't| > 0.74). 
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Figure 10: Measurement of the contamination from uranium noise and electronic noise in the 
e+e~ — > hadrons(7) sample of 1994: The figure shows the transverse energy imbalance for 
event candidates with most of the energy observed either only in the electromagnetic or hadron 
calorimeter and with little matching between tracks in the central tracking chamber and the 
energy deposits. yp. 
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Figure 11: The distribution of the maximum momentum measured in the muon chambers, p^, 
normahzed to the beam energy for e+e~ — > (7) candidates collected in 1993 — 94. The 
contribution from cosmic ray muons is determined from our data as described in the text. 
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Figure 12: The distribution of the acolhnearity angle for e+e — > /x"'"/^ (7) candidates collected 
in 1993 - 94. 
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Figure 13: Rejection of cosmic ray muons in the 1995 data: Figure a) shows the distribution 
of the distance of closest approach to the beam axis, d±. Figure b) shows the absolute value of 
the time, \t\, closest to the beam crossing as measured by the scintillation counters associated 
with the muon candidates. 
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Figure 14: Distribution of the polar angle, defined by the negative muon, of muon pairs collected 
in 1995. The pre-scan and peak data sets are combined. 
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Figure 15: The measured differential cross section e+e~ fi'^fi~{'y) combining the 1993 — 95 
data into three centre-of-mass energy points. The hues show the resuh of a fit using the 
functional form of Equation El 
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Figure 16: Energy of the most energetic jet in e'*"e — > r+r (7) event candidates for 1994 data. 
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Figure 17: Number of tracks associated to each of the two jets in e+e —>■ t^t (7) candidate 
events. 
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Figure 18: The distribution of the acoUinearity for e+e — > t'^t (7) candidates. 
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Figure 19: Distribution of the polar angle of the event thrust axis for e+e" — > t+t'(7) candi- 
dates collected in 1994. The structure seen for | cos^d > 0.65 reflects the modifications of the 
event selection in the end-cap and the transition region. 
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Figure 20: Same as Figure [T31 for e+e r+r (7). 



80 



10 \r eV ^ (y) 
T data 1994 



10 



€"^6 1^1 (y) 




U 1 

I . 



1.1 



E^/E, 



Figure 21: Distribution of the energy of the highest energy cluster measured in the electro- 
magnetic calorimeter for e+e" — > e''"e~(7) candidate events (1994 data). Events below the cut 
value, indicated by the vertical arrow, can be selected by other criteria as described in the text. 
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Figure 22: Acollinearity angle between the two electron candidates in e+e — > e"^e (7) events. 
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Figure 23: The measured differential cross section e+e~ — > e"'"e~(7) for data collected between 
1993 and 1995. The cross sections are calculated for an acoUinearity angle ^ < 25° and a 
mimimum energy of 1 GeV of each final state fermion. The data are compared to the SM 
predictions which are shown as lines for three different centre-of-mass energies. 
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Figure 24: The measured cross sections e+e~ — > hadrons(7) as function of the centre-of-mass 
energy. The sohd hne shows the result of the fit. At the bottom the ratio of the measured cross 
sections and the fit result for the data collected in 1993 — 95 is shown. The errors are statistical 
only. 
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Figure 25: Same as Figure EH for e"'"e (7). 
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Figure 26: Same as Figure 1^ for e"'"e — >• r+r (7). 
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Figure 27: Same as Figure |21] for e+e — >■ e"'"e (7) inside the fiducial volume 44° < 6 < 136° 
for acoUinearity angles ^ < 25° and a minimum energy of 1 GeV of the final state fermions. 
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Figure 28: The measured forward-backward asymmetry in e'^e" — > ix^^~{'^) as function of the 
centre-of-mass energy. The sohd hne shows the result of the fit. At the bottom the difference 
of the measured asymmetry and the fit result for the data collected in 1993 — 95 is shown. The 
errors are statistical only. 
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Figure 29: Same as Figure 1^ for e"'"e — >• r+r (7). 
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Figure 30: Same as Figure 1^ for e^e — > e"'"e (7) in the fiducial volume 44° < 6 < 136°. The 
same cuts as for the total cross section are applied. 
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Figure 31: Distribution of the absolute difference of the measured cross section and forward- 
backward asymmetries (1990 — 95 data) and the five parameter fit result f Table divided by 
the statistical errors of the measurements. The histogram shows the expectation for a Gaussian 
distribution of the measurement in the absence of systematic errors. 
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Figure 32: Contours in the Ap^ — Ri plane for electrons, muons and taus obtained from a fit 
to total cross sections and forward-backward asymmetries. The dashed line shows the contour 
assuming lepton universality and the star indicates the central value. The arrows show the 
change in the SM prediction when varying the input parameters as and mn in the ranges 
defined in Equation ^1 The uncertainty on the SM prediction due to the other parameters is 
small. 
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Figure 33: Contours in the gA-Qw plane for electrons, muons and taus obtained from a fit to 
total cross sections, forward-backward and tau polarisation asymmetries. The dashed line shows 
the contour assuming lepton universality and the star indicates the central value. The arrow 
shows the change in the SM prediction when varying mn in the range defined in Equation ^1 
The uncertainty on the SM prediction due to the other parameters is small. The hollow cross 
indicates the SM expectation without weak radiative corrections. 
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Figure 34: Contour in the JhaV'^z plane obtained from the S-Matrix fit assuming lepton uni- 
versahty. The horizontal band shows the SM predition for with its uncertainty. 
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Figure 35: Contour in the mt-m-w plane obtained from the SM fit to our data. The fit result 
is compared to our direct measurement of the W mass indicated by the hatched band and the 
SM expectation. 
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Figure 36: Contours in the mt-mn plane obtain from the SM fit. The hatched area indicates 
values of the Higgs mass excluded by our direct search result [^2] • 
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Figure 37: The dependence of the SM fit as function of the Higgs boson mass. The shaded 
area indicates the mass range excluded by the direct search. 
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